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ABSTRACT 


The  prediction  of  the  radiated  acoustic  signature  of  underwater  structures  has 
been  and  will  continue  to  be  of  great  importance  for  the  Navy.  New  and  innova¬ 
tive  designs  of  maneuvering  and  propulsion  systems  present  challenging  prob¬ 
lems  for  the  prediction  of  radiated  acoustics.  The  present  research  focused  on  the 
development  of  structural  acoustic  experiments  that  use  air  as  the  working  fluid. 
The  test  specimens  were  designed  using  very  light  weight  structures  such  that  the 
fluid  loading  effects  would  be  non-negligible,  and  in  some  cases  dominant.  A 
number  of  configurations  have  been  investigated,  including  flat  membranes,  rib- 
stiffened  cylindrical  structures,  and  a  ducted  rotor  with  elastic  shell.  The  ducted 
rotor  facility  is  part  of  a  continuing  effort  involving  the  effects  of  inlet  distortion 
on  the  elastic  shell  vibrations.  A  substantial  database  of  flow  field,  acoustic,  and 
vibration  measurements  has  resulted.  The  data  show  a  significant  sensitivity  of 
the  structural  vibration  to  the  spatial  distribution  of  the  acoustic  sources.  In  some 
cases  the  acoustic-vibration  coupling  led  to  a  significant  change  in  the  natural 
frequencies,  and  hence  the  coupling  to  the  frequency  dependent  acoustic  and 
convecting  hydrodynamic  pressure. 
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1.0  INTRODUCTION 

The  ability  to  predict  the  radiated  acoustic  signature  of  underwater  structures  has  been  and  will 
continue  to  be  of  great  importance  for  the  Navy.  New  and  innovative  designs  of  maneuvering  and 
propulsion  systems  present  challenging  problems  for  the  prediction  of  radiated  acoustics  using 
analytical  or  numerical  approaches  alone.  A  significant  obstacle  in  the  development  and  valida¬ 
tion  of  predictive  methodologies  is  a  general  lack  of  high  quality,  comprehensive  experimental 
data.  These  data  are  important  for  both  validation  of  specific  models  and  methodologies,  as  well 
as  for  investigating  phenomena  that  are  not  yet  well  understood. 

The  objective  of  the  present  research  was  to  develop  the  experimental  methods  needed  to  study 
the  coupling  between  hydrodynamics,  hydroacoustics,  and  structural  acoustics  in  an  air-analog 
system.  The  use  of  air  as  a  working  fluid  is  preferred  to  water  in  many  circumstances  due  to  the 
relative  ease  of  working  with  models  and  instrumentation  in  a  laboratory.  An  air  environment 
allows  for  standard  laboratory  materials  and  instrumentation  to  be  used  without  difficulty.  Addi¬ 
tionally,  models  in  air  can  be  placed  in  a  near-perfect  anechoic  environment  in  order  to  easily 
obtain  radiated  sound  power,  directivity,  etc. 

The  relevant  parameters  for  fluid-flow  problems  with  rigid  boundary  conditions  include  the  Rey¬ 
nolds  number  and  the  Mach  number.  This  allows  most  hydrodynamic  and  hydroacoustic  problems 
to  be  considered  in  air  with  little  difficulty.  The  use  of  air  for  the  study  of  underwater  structural 
vibrations  and  acoustics  is  generally  not  considered  acceptable  since  dimensional  similarity  is 
generally  not  possible.  Since  the  materials  that  can  typically  be  used  for  air  experiments  are  not 
too  different  from  those  used  in  underwater  applications,  there  is  a  necessary  mismatch  between 
the  dimensionless  parameters  of  an  underwater  application  and  an  air-analog  experiment.  How¬ 
ever.  the  effects  of  membrane  waves,  flexural  waves,  fluid  loading  can  be  studied  in  air-analog 
systems  under  parameter  ranges  that  are  of  direct  relevance  to  underwater  problems.  For  example, 
it  will  be  shown  that  a  thin  plate  or  shell  exhibits  significant  fluid  loading  effects  in  air  at  tow  fre¬ 
quencies.  By  considering  higher  frequencies  or  thicker  materials,  the  effects  of  dispersive  waves 
with  little  fluid  loading  can  be  examined.  This  flexibility  in  designing  the  air-analog  experiments 
presents  a  unique  advantage  over  water  experiments. 

A  number  of  various  geometries  were  considered  in  this  study.  The  intent  was  to  gain  experience 
with  the  methods  and  instrumentation  using  relatively  simple  problems,  and  then  consider  config- 
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urations  of  increasing  complexity.  The  majority  of  the  efforts  were  focused  on  the  development 
and  measurement  of  a  new  facility  that  incorporates  a  rotor  operating  within  an  elastic  duct.  The 
geometry  was  considered  to  have  a  suitable  balance  between  simplicity  in  its  boundary  conditions 
and  complexity  of  the  physical  phenomena  involved.  Other  geometries  included  combinations  of 
plates  and  shells  with  both  hydrodynamic  and  acoustic  excitation. 

It  is  noted  that  the  intended  focus  of  this  research  evolved  over  the  course  of  the  project.  Discus¬ 
sion  w  ith  engineers  from  the  NSWC,  ONR,  and  Electric  Boat  allowed  the  experimental  program 
to  take  advantage  of  a  number  of  capabilities  that  were  of  direct  interest  to  the  Navy.  Specifically, 
the  effects  of  internal  mean  flow  distortion  on  the  structural  response  of  an  elastic  duct  was  con¬ 
sidered  in  some  detail. 

The  remainder  of  this  report  is  structured  as  follows.  Section  2  briefly  describes  the  basic  instru¬ 
mentation  that  was  used  in  this  research.  Section  3  provides  a  detailed  description  of  the  geometry 
and  properties  of  the  test  setup  and  the  elastic  structures  used.  Sections  4  though  8  describe  the 
primary  results  of  the  experimental  program  including  a  detailed  set  of  conclusions  in  each  sec¬ 
tion.  Overall  conclusions  and  a  discussion  of  the  project  are  given  in  Section  9. 
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2.0  EXPERIMENTAL  TECHNIQUES 

The  experiments  described  in  the  following  sections  of  this  report  will  require  a  number  of 
sophisticated  measurement  systems.  A  summary  of  the  measurement  capabilities  and  relevant  ref¬ 
erences  are  given  in  the  following  sub-sections. 

2.1  Anechoic  Wind  Tunnel 

All  experiments  will  be  performed  in  the  Anechoic  Wind  Tunnel  (AWT)  located  within  the  Hes- 
sert  Laboratories  for  Aerospace  Research  at  the  University  of  Notre  Dame.  A  schematic  of  the 
tunnel  is  shown  in  Figure  I .  The  AWT  consists  of  an  anechoic  chamber  w  ith  a  working  space  that 
is  20ft  w  ide  by  26ft  long  by  8ft  high.  The  22 in  fiberglass  wedges  have  a  coefficient  of  energy 
absorption  at  normal  incidence  of  0.99  or  greater.  The  wedges  provide  a  low  frequency  cutoflfof 
about  100Hz. 

A  low  disturbance  subsonic  free-jet  test  section  wind  tunnel  is  located  within  the  anechoic  cham¬ 
ber.  The  cross-sectional  area  of  the  test  region  is  4ft2.  A  low  noise  centrifugal  blower  is  connected 
to  the  flow-exit  of  the  tunnel  (not  shown  in  Figure  1 )  through  a  muffler/diffuser  system.  This  pro¬ 
vides  a  maximum  test  section  velocity  of  1 1  Oft/sec.  A  complete  description  of  the  tunnel's  indi¬ 
vidual  components  and  important  characteristics  is  presented  in  Mueller  et  al.  (1992a.b). 
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2.2  Phased  Microphone  Array 

The  AWT  is  equipped  with  a  dual  planar  array  measurement  system  for  acoustic  measurements. 
Measuring  acoustics  using  an  array  of  microphones  offers  two  major  advantages  over  measuring 
acoustics  using  a  single  microphone.  First,  the  array  increases  the  signal  to  noise  ratio  over  that  of 
a  single  microphone  by  a  factor  approximately  equal  to  the  number  of  microphones.  The  second 
benefit  of  an  array  is  its  capability  to  spatially  localize  acoustic  sources.  The  two,  40  microphone 
arrays  are  fully  calibrated,  and  capable  of  reducing  the  ambient  noise  levels  of  the  facility  by  as 
much  as  20dB. 

F.ach  of  the  AWT  planar  arrays  consists  of  a  strategic  arrangement  of  40  microphones  mounted  in 
a  flat  panel.  Figure  2  shows  a  picture  of  one  array  panel  situated  outside  of  the  AWT  before 
microphone  installation,  and  a  completed  acoustic  array  situated  inside  the  AWT  from  the  back¬ 
side  of  the  array  panel.  By  placing  microphones  to  avoid  replicating  the  spacing  between  any  pair, 
an  acoustic  array  resolves  the  acoustic  contributions  from  a  focal  location  while  rejecting  acoustic 
contributions  from  other  locations  {Underbrink.  2002).  By  using  perforated  plastic  sheets  as  the 
array  panels,  the  locations  of  the  microphones  in  the  AWT  array  panels  may  be  altered  to  suit  the 
needs  of  the  acoustic  measurement  experiments  at  hand.  Also,  all  microphones  from  the  two 
arrays  may  be  merged  into  one  array  panel  for  single  plane  measurements. 

For  the  sake  of  example.  Figure  3  shows  a  diagram  of  the  acoustic  array  panels  in  the  two-panel 
configuration  with  an  array  located  on  either  side  of  the  open  jet  test  section  to  measure  the  trail¬ 
ing  edge  acoustics  from  a  model  foil.  In  the  model  foil  experiments,  the  trailing  edge  is  oriented 
vertically  in  the  wind  tunnel  and  the  open  jet  is  constrained  by  acoustically  absorbing  endplates 
that  face  the  model  side  edges. 


Figure  2  Picture  of  phased  array  panel  (left),  and  completed  inside  the  AWT  (right) 
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In  the  acquisition  stage  of  the  acoustic  array  experiments,  a  digital  representation  of  the  analog 
signal  from  each  microphone  is  acquired  using  10  synchronized  analog  to  digital  converters  con¬ 
trolled  by  a  PC  workstation.  After  time  series  data  have  been  acquired  from  ail  the  microphones, 
post-processing  algorithms  focus  the  array  to  a  grid  of  locations  in  space.  The  most  common  array 
source  localization  technique  is  a  post-processing  operation  called  beam  forming. 

Beamforming  produces  an  acoustic  source  level  map  over  a  specified  grid  of  points  in  space.  The 
classical  beam  forming  process  shifts  each  microphone  signal  in  time  by  a  delay  equal  to  the 
assumed  acoustic  wave  travel  time  from  the  beamforming  grid  point  to  the  microphone  location. 
The  shifted  microphone  signals  are  summed  either  in  the  time  or  the  frequency  domain.  The  mean 
square  value,  or  the  power  spectral  density  level,  represents  a  value  proportional  to  the  mean- 
square  acoustic  pressure  emanating  from  the  beam  forming  grid  point.  Scanning  a  grid  of  points  in 
space  forms  a  beamforming  map.  I  fan  acoustic  source  exists  at  the  beam  forming  grid  point,  the 
summed  acoustic  source  signal  is  amplified  and  a  peak  appears  in  the  beamforming  map.  In  addi¬ 
tion.  extraneous  noise  that  is  uncorrelated  among  microphones  is  reduced  through  phasing  the 


AWT  Test  Section 
Looking  upstream 


End  plate  oource 

Figure  3  Schematic  of  dual  array  configuration  inside  the  AWT 


array  via  beamforming.  The  array  is  calibrated  using  a  speaker  to  convert  the  beamforming  map 
into  a  sound  pressure  level  (SPL)  map.  Figure  4  shows  an  SPl.  map  produced  via  beam  forming 
using  measurements  in  the  presence  of  a  calibration  speaker  for  a  frequency  of  3kHz. 


A  full  description  of  classical  beamforming  in  the  frequency  domain  follow's  the  work  of  Dough¬ 


erty  (2002).  Beamforming  has  two  limitations,  which  are  minimized  through  array  pattern  design 
optimization.  The  first  is  the  capability  of  the  array  to  resolve  independent  sources.  The  array  res¬ 
olution  is  quantified  by  the  array  main  lobe  (or  beam),  which  depends  on  the  array  aperture,  the 
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outer  perimeter  of  the  array,  relative  to  the  wavelength  of  the  sound  and  the  distance  from  the 
sound  source  to  the  array.  Figure  4  points  out  the  main  lobe  in  the  beamforming  map  of  the  cali¬ 
bration  speaker  plane.  The  larger  the  wavelength  and  the  larger  the  distance  from  the  sound 
source  to  the  array,  the  wider  the  beam  width.  The  larger  the  array  aperture,  the  more  narrow  the 
beam  width.  The  number  of  microphones  and  their  placement  governs  the  second  limitation  of  the 
array:  the  strength  of  sidelobes.  Sidelobes  are  unwanted  artifacts  of  beamforming  that  appear  as 
secondary  peaks  in  the  beam  forming  map  as  illustrated  in  Figure  4.  They  are  essentially  the  prod¬ 
uct  of  spatial  ambiguity  due  to  2 n  multiples  in  phase  among  more  widely  spaced  microphone 
pairs  in  the  presence  of  a  source.  The  nearness  of  the  sidelobes  to  the  mainlobe  and  the  magnitude 
of  the  sidelobes  increase  with  decreasing  wavelength  and  increasing  closeness  to  the  source. 


BMmtamllv,  3203HI.  SPL  <<S  raf  IlfcFa) 


Figure  4  S  PL  map  produced  via  beam  forming  using  measurements  from  the  array  in  AWT 

2.3  Laser  Doppler  Vibrometry 

A  laser  Doppler  vibrometer  is  based  on  the  principle  of  the  detection  of  the  Doppler  shift  of 
coherent  laser  light  that  is  scattered  from  a  small  area  of  the  test  object.  The  object  scatters  or 
reflects  light  from  the  laser  beam  and  the  Doppler  frequency  shift  is  used  to  measure  the  compo¬ 
nent  of  velocity  which  lies  along  the  axis  of  the  laser  beam.  Since  the  laser  light  has  a  very  high 

frequency  (typically  >  1 0 1  ^Hz)  a  direct  demodulation  of  the  light  is  not  possible.  An  optical  inter¬ 
ferometer  is  therefore  used  to  mix  the  scattered  light  coherently  with  a  reference  bream.  The 
photo  detector  measures  the  intensity  of  the  mixed  light  whose  beat  frequency  is  equal  to  the  dif¬ 
ference  frequency  between  the  reference  and  the  measurement  beam.  Vibration  frequencies  up  to 
400kHz  can  be  measured  from  a  distance  of  up  to  20m. 
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Two  LDV  systems  were  used  for  the  structural  vibration  measurements.  The  first  is  a  Polytec  sin¬ 
gle  point  LDV  system.  The  second  is  a  Poly  tec  Scanning  system  that  allows  for  fully  program  able 
geometry  scans  for  a  given  experiment.  The  two  systems  were  used  in  combination  in  order  to 
obtain  the  full  two-point  correlation  functions  of  vibrational  velocity.  This  information  is  impor¬ 
tant  when  stochastic  forcing  such  as  when  broadband  rotor  noise  and/or  the  convening  hydrody¬ 
namic  pressure  is  dominant. 
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3.0  EXPERIMENTAL  APPARATUS 


3.1  Tensioned  Plate/Membrane  Apparatus 

A  new  experiment  was  created  in  which  thin  materials  could  be  placed  in  the  AWT  with  arbitrary 
tension.  The  purpose  was  to  use  acoustic  and  hydrodynamic  forcing  in  order  to  gain  experience 
working  with  the  structural  acoustics  in  very  thin  materials.  This  section  will  describe  the  geome¬ 
try  and  the  acoustic  response  experiments  using  a  very  thin  aluminum  membrane.  Section  6  of 
this  report  will  describe  the  response  of  the  membrane  to  complex  hydrodynamic  forcing. 

A  schematic  of  the  fixture  geometry  is  shown  in  Figure  5.  The  assembly  consists  of  an  aluminum 
frame  within  a  larger  structure.  The  inner  frame  supports  the  test  specimen.  An  insert  piece  is  used 
to  define  the  shape  of  the  vibrationally  active  region  of  the  membrane.  The  material  is  tensioned 
from  four  sides  using  a  mechanical  pivot  with  a  spring  and  force  transducer.  Front  and  back 
images  of  the  completed  assembly  in  the  AWT  are  shown  in  Figure  6.  The  active  region  was  cir¬ 
cular  as  observed  on  the  image  of  the  back  of  the  fixture.  The  material  used  in  the  data  presented 
in  this  section  was  aluminum  with  a  thickness  of  0.0127mm  (0.0005  inches).  The  mass  of  the 
membrane  per  unit  surface  area  was  equivalent  to  2.8cm  of  air. 

The  intent  of  the  exercise  was  to  measure  the  natural  frequencies  and  mode  shapes  of  the  mem¬ 
brane  using  various  materials  under  various  tension  values.  For  a  circular  membrane  the  motion 
normal  to  the  plane  is  a  solution  to 


LAf  +  l  5J1  =  1  a2n 

rdr\  dr)  c 2  dt2 


0) 


where  c2  =  T/(p$h)  is  the  wave  speed  of  the  membrane  as  set  by  the  Tension  (T)  and  the  mass 


per  unit  area  (p  sh).  The  discrete  solutions  of  equation  1  are  classified  by  the  mode  order  pair  (m.n) 
where  m  is  the  number  of  diametrical  nodal  lines  and  n  is  the  number  of  circular  nodal  lines 
including  the  outer  boundary  of  the  material.  The  first  several  mode  shapes  are  shown  in  Figure  7. 
The  vibration  of  the  membrane  was  measured  using  the  single  SLDV.  The  measurement  grid 
computed  as  shown  in  Figure  8.  These  points  were  programed  into  the  SLDV  lor  automatic  scan¬ 
ning.  The  membrane  was  forced  using  broadband  acoustic  forcing  and  a  small  speaker  placed  in 
proximity  to  the  structure.  Two  examples  that  identify  the  (1,1 )  and  (0,2)  modes  are  shown  in  Fig¬ 
ure  9.  The  colors  are  representative  of  the  amplitude  of  the  motion  at  a  particular  phase. 
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The  autospectral  density  from  each  of  the  measurement  points  were  averaged  in  order  to  provide 
the  net  vibration  of  the  membrane.  These  data,  along  with  a  number  of  identified  mode  order  pairs 
are  shown  in  Figure  10  using  a  tension  of  1 82.5N/m.  The  natural  frequencies  can  be  easily 
observed  by  the  local  maxima  in  this  spectrum.  The  full  two  point  correlations  were  evaluated  at 
each  of  these  maximum  (e.g.,  figure  9)  in  order  to  identify  the  mode  order  pair.  Table  I  summa¬ 
rizes  the  measured  and  calculated  frequencies  of  the  first  few  modes.  It  can  be  observed  that  the 
measured  natural  frequencies  are  lower  than  the  in-vacuo  values  due  to  the  significant  fluid  load¬ 
ing  effects.  For  example,  the  lowest  mode  (0, 1 )  indicates  a  resonant  frequency  of  36 Hz,  compared 
to  the  in-vacuo  value  of  137Hz.  The  reduction  in  frequency  by  a  factor  of  3.8  indicates  that  the 
acceleration  of  the  mass  of  the  air  surrounding  the  membrane  is  approximately  14.5  times  the 
mass  of  the  membrane  itself  (i.e.,  the  square  of  the  frequency  ratio).  The  effective  thickness  of  the 
accelerated  air  on  each  side  of  the  membrane  is  also  given  in  Table  I . 

It  is  of  interest  to  note  that  the  fluid  loading  factor  (Blake,  1986).  defined  as: 


(pc) 


'Jluhl 


(f>sh)  a 

is  approximately  50  at  the  lowest  mode  (36Hz).  That  is.  despite  the  relatively  large  density  ratio 


(2) 


from  the  aluminum  to  air  (2700/1.15),  the  structural  vibrations  are  dominated  by  the  fluid  loading 
effects  because  the  product  of  the  acoustic  wavenumber  and  the  material  thickness  is  also  large.  In 
this  example  (36 Hz),  knh=  \  1 8,000. 


The  ratio  of  the  measured  natural  frequency  of  each  mode  to  the  in-vacuo  natural  frequency  can 
be  observed  to  increase  for  each  higher  order  mode.  This  is  expected  since  the  fluid  loading  factor 
is  less  for  each  successive  mode.  Results  from  a  larger  tension  value  of255.4N/m  are  shown  in 
Figure  1 1  and  Table  2.  These  data  show  similar  results  in  terms  of  the  changes  in  natural  fre¬ 
quency.  The  frequency  ratio  { fmea/fin-vacuo ) 's  plotted  in  Figure  12  for  both  sets  of  data.  The  fre¬ 
quencies  have  been  further  normalized  by  the  ratio  of  the  membrane  wave  speed  to  the  sound 
speed  (c„/c0). 
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Figure  5  Schematic  of  mechanical  fixture  for  plate/membrane  stretching 


Figure  7  Mode  shapes  of  a  circular  membrane  (from  Morse  and  Ingard) 
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Figure  8  Scanning  points  on  circular  active  region  of  aluminum  membrane 


Figure  9  Example  of  measured  eigenmode  shape  of  circular  membrane  for  the  (1,1)  mode  at 
78Hz  (left)  and  the  (0,2)  mode  at  1 85Hz  (right). 
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10  Spatially  averaged  auto  spectral  density  of  the  membrane  vibration  (T=18 

Frequency 

Mode 

In  Vacuo 

Frequency 

Amount  of 

(Hz) 

(tfvi) 

Frequency  (Hz) 

Ratio 

Air  (cm) 

36 

(0,1) 

137.38 

0.26204 

20.381 

78 

(U) 

218.89 

0.35633 

11.022 

127 

(2.1) 

293.28 

0.43175 

7.5275 

170 

(3,1) 

364.48 

0.46641 

6.4335 

185 

(0,2) 

315.35 

0.58665 

4.0665 

241 

(0,3) 

494.36 

0.48750 

5.8892 

284 

(2,2) 

480.85 

0.59062 

4.0147 

Table  i:  Frequencies  and  corresponding  mode  numbers  for  fluid  loaded  foil  at  182.5N/m  com¬ 
pared  to  the  in  vacuo  frequencies. 
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Figure  1 1  Spatially  averaged  auto  spectral  density  of  the  membrane  vibration  (T=255.4N/m) 


Frequency  Mode  In  Vacuo  Frequency  Amount  of 

(Hz) (jw,fi)  Frequency  (Hz) Ratio_ Air  (cm) 


38 

(0,1) 

162.55 

0.23377 

25.610 

83 

(U) 

259.00 

0.32046 

13.628 

128 

(0,2) 

373.12 

0.34305 

11.893 

190 

(1,2) 

474.21 

0.40067 

8.7184 

244 

(0,3) 

584.94 

0.41714 

8.0533 

245 

(2,2) 

568.95 

0.43062 

7.5478 

306 

0,3) 

687.67 

0.44426 

7.0983 

358 

(0.4) 

797.03 

0.45607 

6.7528 

Table  2:  Frequencies  and  corresponding  mode  numbers  for  Hu  id  loaded  foil  at  255,4N/m  com¬ 


pared  to  the  in  vacuo  frequencies. 
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Figure  12  Ratio  of  the  observed  natural  frequency  to  the  in-vacuo  frequency,  normalized  by 
the  wave  speed  ratio  ( cn/ca )  as  a  function  of  the  fluid  loading  factor  (p).  The  solid  line  repre¬ 
sents  the  function  0.27(P‘a4)  is  shown  only  for  visual  reference. 
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3.2  Nickel  Shell  Apparatus 

One  of  the  significant  results  from  this  project  is  the  development  of  a  set  of  cylindrical  test  spec¬ 
imens  that  are  very  thin,  yet  do  not  have  welds,  seams,  or  other  large  scale  imperfections  that 
result  from  the  manufacturing  process.  The  purpose  in  developing  these  test  shells  was  to  provide 
the  symmetry  needed  for  simplified  analysis  of  the  structure. 

The  manufacture  of  the  thin  shells  was  as  follows.  First,  an  acrylic  cylinder  with  8. 1  inch  outer 
diameter  was  machined.  The  elastic  shell  was  then  created  by  electroplating  pure  nickel  onto  the 
acrylic.  Note  that  a  minute  layer  of  gold  was  first  sputtered  onto  the  base  cylinder  in  order  to  pro¬ 
vide  the  electrical  conductivity  needed  for  the  electroplating  process.  The  time/current-density 
were  monitored  to  provide  a  material  that  is  essentially  flawless  and  0.0027  inches  in  thickness. 
Support  ring  (ribs)  were  mounted  and  glued  to  the  nickel  shell  for  structural  support.  Axial  sup¬ 
port  struts  were  then  placed  betw-een  the  ribs  in  order  to  prevent  yielding  or  buckling  of  the  thin 
shell  during  handling.  The  final  step  of  the  process  required  the  submersion  of  the  acrylic/nickel 
shell  into  an  acidic  solution  in  order  to  chemically  remove  the  acrylic  base  material.  Once  all  of 
the  inner  material  was  removed,  the  final  product  was  removed  from  the  solution  and  cleaned. 

A  total  of  four  nickel  duct  sections  were  constructed.  Two  images  of  one  of  these  sections  is 
shown  in  Figure  13.  The  supporting  apparatus  shown  provides  the  capability  to  pre-stress  the 
material  with  axial  tension  without  applying  unwanted  torsional  pre-stress.  The  system  uses  linear 
bearings  and  two  ground  rods  that  support  the  rib  sections.  Four  linear  springs  are  compressed  in 
order  to  provide  a  desired  axial  tension.  The  linear  bearings  attached  to  the  spring  section  are  held 
in  place  (axially)  using  two  force  transducer  buttons  that  provide  an  accurate  measurement  of  the 
applied  tension. 
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33  Dueled  Rotor  System 

A  new  ducted  rotor  facility  was  developed  for  the  experiments  that  will  be  described  in  later  Sec¬ 
tions.  The  intent  of  the  facility  was  to  create  boundary  conditions  for  the  flow-structu  re-acoustic 
problem  that  captured  the  complexity  of  the  physical  mechanisms  involved  with  these  types  of 
engineering  systems,  while  maintaining  the  symmetry  and  simplicity  that  allows  the  problem  to 
be  studied  analytically. 

A  schematic  of  the  ducted  rotor  system  with  rigid  walls  is  shown  in  Figure  14.  The  duct  was  con¬ 
structed  from  PVC  with  a  machined  inner  diameter  of  206  mm  and  a  wall  thickness  of  6.7  mm. 
and  supported  by  a  vibration  isolation  pad.  A  center  cylinder  extended  upstream  of  the  duct  inlet 
and  downstream  of  the  outlet  to  create  an  annular  flow  passage.  The  inlet  flow  was  conditioned 
using  a  single  layer  of  a  "cheesecloth"  type  thin  fabric  with  negligible  acoustic  impedance 
stretched  over  a  cylindrical  wire  frame.  A  compact  servo  motor  was  housed  in  the  center  cylinder 
to  direct  drive  the  rotor  at  the  desired  speed.  The  upstream  cylinder  and  motor  were  supported  by 
eight  streamlined  struts.  The  downstream  center  cylinder  was  cantilevered  from  a  support  struc¬ 
ture  beyond  the  duct  outlet  (not  shown).  The  rotor  was  a  ten -b laded  propeller  previously  used  by 
Sevik,  with  diameter  D  =  0.203  m,  constant  blade  chord  C  =  0.025  m  and  maximum  thickness  of 
2.4  mm.  with  straight  blades  incorporating  twist  only.  The  rotor  operated  with  a  tip  gap  equal  to 
5%  of  the  blade  chord.  The  length  of  the  duct  was  variable,  as  was  the  position  of  the  rotor  from 
the  duct  inlet. 

The  turbulent  boundary  layer  created  on  the  sidewalls  of  the  duct  upstream  of  the  rotor  plays  an 
important  role  in  the  modeling  and  prediction  of  the  acoustic  sources.  The  velocity  field  was  char¬ 
acterized  using  a  single  component  hot-wire  anemometer  connected  to  a  radial  traverse  system. 
Time  series  of  the  velocity  were  measured  at  a  range  of  radial  and  axial  positions.  The  mean  and 
RMS  of  velocity  from  these  measurements  is  shown  in  Figure  15.  The  data  indicate  the  expected 
acceleration  effects  in  the  very-near  entrance  region,  followed  by  a  typical  boundary  layer  that 
grows  smoothly  in  the  downstream  direction.  These  data  will  be  used  in  later  sections  in  defining 
the  acoustic  source  strength  with  both  rigid  and  elastic  shells. 

Example  sound  spectra  are  shown  in  Figures  16  and  17.  A  number  of  distinctive  features  were 
observed  in  the  spectra.  For  example,  multiple  "organ  pipe"  type  resonances  can  be  observed  at 
fixed  frequencies  in  Figure  16.  The  spectra  are  plotted  with  frequency  normalized  by  the  blade 
passing  frequency  in  Figure  17.  These  data  show  distinct  features  at  integer  multiples  of  the  blade 
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rate.  Significant  broad  band  sound  was  also  observed  at  high  frequencies.  The  radiated  sound  was 
found  to  be  dependent  on  the  rotor  flow  coefficient.  The  duct  flow  rate  was  decreased  by  applying 


variable  restriction  in  the  form  of  thin  fabric  with  negligible  effect  on  the  acoustic  field  to  the  duct 


exit,  and  increased  by  attaching  the  outlet  of  the  duct  to  the  anechoic  wind  tunnel  blower. 


The  rotor  casing  and  duct  were  segmented  such  that  individual  sections  could  be  removed  and 
replaced  with  one  or  more  of  the  elastic  shell  sections.  Model  schematics  of  the  ducted  rotor  with 
the  nickel  shell  and  tensioning  apparatus  are  shown  in  Figure  18.  The  right-most  picture  shows 
the  shell  with  a  rigid  walled  exit  contraction.  This  configuration  allows  the  internal  pressure  of  the 
duct  to  be  controlled.  Specifically,  the  flow  rate  generated  by  the  rotor  is  accelerated  through  the 
contraction  and  exits  as  a  jet  at  atmospheric  pressure.  The  acceleration  therefore  provides  an  inter¬ 
nal  mean  pressure  that  can  be  controlled  by  either  changing  the  flow  rate  through  the  duct,  or  by 
changing  the  contraction  ratio.  This  feature,  along  with  the  variable  tension  provided  by  the  sup¬ 
porting  mechanism,  allows  for  some  control  of  the  shell's  natural  frequencies  by  changing  the 
static  pre-stress. 
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Support  Struts 


Figure  14  Schematic  of  the  ducted  rotor  system 
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Figure  15  Characterization  of  the  duct  inflow  boundary  layer  by  the  mean  velocity  pro¬ 
files  (left)  and  the  RMS  of  velocity  profiles  (right).  The  location  numbers  noted  in  the 
legend  represent  the  distance  (in  inches)  from  the  inlet. 


Frequency,  Hz 

Figure  1 6  Measured  sound  spectra  for  the  ducted  rotor  at  three  rotational  speeds. 
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Figure  1 7  Measured  sound  spectra  for  the  ducted  rotor  with  the  frequency  normalized  by  the 

blade  passing  frequency. 


Figure  18  Model  schematic  of  ducted  rotor  with  nickel  section  and  tensioning  device  (left),  and 

rotor/she  1 1  with  exit  contraction  (right). 
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4.  RECIPROCITY  MEASUREMENT  OF  ACOUSTIC  TRANSFER  FUNC¬ 
TIONS 

The  prediction  of  structurally  and  aerodynamically  generated  sound  requires  detailed  knowledge 
of  both  the  distribution  of  sources  as  well  as  the  transfer  of  acoustic  energy  from  the  source  to  a 
location  of  interest.  The  propagation  of  acoustic  waves  can  be  found  from  analytical  methods  or 
numerical  solutions  when  the  boundary  conditions  are  sufficiently  simple.  However,  complex 
geometries  involved  with  engineering  systems  can  result  in  significant  difficulty  for  these 
approaches,  particularly  when  surfaces  are  vibrationally  active,  when  sound  absorbing  materials 
are  used,  or  when  features  of  the  geometry  are  complicated.  The  present  communication  describes 
a  reciprocal  experiment  using  a  miniature  dual -microphone  probe.  The  measurements  allow  the 
calculation  of  the  frequency  dependent  Green's  function  and  its  spatial  gradient,  which  represent 
the  acoustic  transfer  functions  for  monopole  and  dipole  sources,  respectively. 

A  point  monopole  with  a  given  frequency  (to)  can  be  represented  by  its  complex  source  strength, 
5(.rs.) .  This  results  in  an  acoustic  pressure  at  location  x„  that  is  given  by  the  solution  to  the  I  lei- 
moltz  equation  (Pierce,  1989): 

(V2  +  k2)p(x0)  =  -4ic5(x()5(x0  -*.) .  (3) 

where  k  is  the  acoustic  wavenumber.  This  expression  is  valid  if  the  fluid  is  at  rest,  and  approxi¬ 
mately  true  at  low  Mach  number.  The  Green’s  function  G(x0,  x  )  is  defined  as  the  solution  to 
equation  (3)  for  a  unit  source  strength  (S=l).  The  acoustic  pressure  can  then  be  determined  by: 

p(x0)  =  G(x0,xs)S(xs).  (4) 

That  is,  the  function  G  represents  the  “transfer  function"  between  the  acoustic  source  at  location 
xs  and  the  acoustic  pressure  at  location  x0. 

Similarly,  consider  a  single  dipole  with  vector  valued  amplitude  .  The  inhomogeneous 
Helmholtz  equation  can  be  written  as 

(V2  +  k2)p(x0 )  =  -4jiZ)  ■  V8(x-xs) .  (5) 

The  acoustic  pressure  field  is  given  as 

P(x0)  =  t>(xs)  ■  VG(x0,xs) .  (6) 
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1 1  will  be  useful  to  consider  a  dipole  with  an  assumed  fixed  vector  direction,  such  that  a  new  func¬ 
tion  can  be  defined  as: 

Gl}(xo,xs)^0VG)/\ti\  .  (7) 

The  radiated  sound  can  now  be  given  by 

P(xa)  =  Gd(xo>xs)  |£(-v,)|  (8) 

The  function  G^xvxs )  therefore  represents  the  acoustic  transfer  function  for  dipoles  with  a  fixed 
direction.  This  representation  is  particularly  useful  in  applications  where  aeroacoustic  source  dis¬ 
tributions  are  the  subject  of  interest.  The  specific  example  to  be  considered  later  in  this  paper  is 
that  of  a  short,  three  dimensional  duct  section.  A  low-speed  rotor  (such  as  a  fan  or  compressor) 
operating  in  such  a  duct  section  will  create  a  distribution  of  dipole- 1  ike  sources  that  lead  to  far 

field  sound.  Estimates  of  |z)|  can  be  obtained  in  this  example  by  calculating  the  unsteady  lift  on 
the  blades  due  to  approach  turbulence  (see,  e.g.,  Goldstein.  1976).  Equation  6  (or  its  integral  when 

Z$  is  represented  as  a  force  per  unit  volume)  provides  a  simple  method  for  predicting  the  far  field 
sound  pressure  level  when  GD  is  known. 

The  duct  and  rotor  combination  described  can  also  be  considered  as  an  acoustic  inverse  problem. 
Specifically,  acoustic  measurements  can  be  made  using  a  single  microphone  located  in  the  far 
field.  In  this  example,  an  independent  estimate  of  the  acoustic  transfer  function,  G'^,  will  allow 
the  determination  of  the  net  dipole  source  strength. 

A  direct  measurement  of  G  would  require  the  placement  of  a  monopole  type  source  at  a  desired 
location  (.vv),  and  a  microphone  located  at  the  desired  observation  location  (*„).  In  practice,  how¬ 
ever,  it  is  not  generally  possible  to  create  an  omnidirectional  point  monopole,  particularly  w'hen 
the  geometry  of  interest  is  small  or  difficult  to  access.  Furthermore,  the  source  strength  of  such  a 
device  could  be  modified  by  strong  near  field  pressures  resulting  from  rigid  and  vibrating  surfaces 
in  the  region  of  the  source.  The  direct  measurement  of  GD  would  require  a  known  amplitude  ideal 
dipole  to  be  placed  at  xs  which  is  not  feasible  in  practice. 

A  method  for  circumventing  these  issues  is  to  use  a  reciprocal  technique.  The  principal  of  acous¬ 
tic  reciprocity  was  first  proposed  by  Helmholtz  in  I860,  and  later  given  by  Rayleigh  in  1873.  A 
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general  proof  of  these  concepts  was  later  given  by  Lyamshev  (1959)  which  utilized  Green's  for¬ 
mula  to  show  that  reciprocity  was  valid  for  linear,  self-adjoint  differential  operators.  Since  that 
time,  a  number  of  authors  have  provided  further  refinements  to  both  the  theory  and  experimental 
application.  For  example,  Steenhoek  and  Wolde  (1970)  and  Belousov  and  Rimskii-Korsakov 
(1975  )  both  provide  detailed  discussions  on  the  application  of  reciprocity  to  mechanical-acousti¬ 
cal  systems.  Fahy  (1995.  2003)  provided  a  review  of  the  development  of  reciprocity  as  well  as 
detailed  schematics  of  the  reciprocal  relations  for  a  large  number  of  vibroacoustic  scenarios. 

The  use  of  reciprocity  for  aeroacoustic  applications  is  significantly  more  rare.  Presumably  this  is 
because  of  the  multi-pole  nature  of  rigid  body  aeroacoustic  problems.  In  the  papers  by  Fahy 
( 1995,  2003)  the  reciprocity  equations  for  both  monopole  and  dipole  point  sources  embedded  in  a 
fluid  were  given.  It  is  these  relations  that  will  allow  for  the  experimental  determination  of  the 
functions  O’  and  Gq.  To  the  author's  knowledge,  this  communication  presents  the  first  application 
of  Fahy's  dipole  form  of  acoustic  reciprocity  for  experimental  measurements. 

Although  the  basic  acoustic  relationships  described  in  this  paper  are  well  understood,  the  experi¬ 
mental  methodology  and  associated  uncertainty  considerations  are  not.  The  techniques  described 
herein  show  that  relatively  simple  experiments  can  be  conducted  to  determine  both  G  and  GD 
simultaneously.  The  measurement  uncertainty  will  be  first  considered  analytically,  and  then  dem¬ 
onstrated  using  results  obtained  in  a  free  space  (anechoic)  environment.  An  example  experiment 
will  also  be  described  which  illustrates  the  use  of  the  technique  in  measuring  the  monopole  and 
dipole  acoustic  propagation  from  inside  of  a  short,  rigid  duct  section  to  a  location  outside  the  duct. 

4.1  Application  of  monopole  and  dipole  reciprocity 

A  schematic  of  monopole  and  dipole  reciprocity  adapted  from  Fahy  (1995)  is  shown  in  Figure  19. 
The  geometry  is  arbitrary,  and  shown  here  as  a  short  cylindrical  duct  as  it  pertains  to  the  example 
used  later  in  this  paper.  For  point  monopoles,  the  spatial  locations  of  the  source  and  the  micro¬ 
phone  can  simply  be  switched.  The  Green's  function  as  defined  in  equation  4  can  then  be  calcu¬ 
lated  as: 

P(Xo>  m  Pits) 

S(xs)  S(x0) 
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where  the  term  on  the  left  represents  the  desired  Green's  function,  and  the  term  on  the  right  is  the 
reciprocal  representation  that  is  measured  by  placing  a  microphone  at*s  and  a  calibrated  omnidi¬ 
rectional  acoustic  source  at  jc0. 

The  reciprocal  relationship  for  dipoles  is  given  as: 


»(-Ol  =  |fl(X,)| 
t)(xs)  |25(*0)|  ‘ 


(10) 


That  is,  the  ratio  of  the  acoustic  particle  velocity  magnitude,  |/I(jc0)| ,  to  the  dipole  source  magni¬ 
tude.  p(xf)| ,  is  fixed  when  the  locations  .v„  and.r,  are  switched.  This  is  true  only  if  the  directions 
of  the  velocity  and  dipole  vectors  are  maintained  as  shown  in  Figure  19. 

Two  important  simplifications  of  equation  8  can  be  made  if*,,  is  located  in  the  far  field.  First,  the 
acoustic  velocity  can  be  represented  by 

«(x0)  =  p(xa)/pc.  (11) 

Second,  for  wavenumber  vectors  parallel  to  the  line  segment  connecting x0  and  xs,  the  sound  field 
incident  upon  the  geometry  from  a  dipole  located  atjc„  is  identical  to  that  of  a  monopole  whose 
strength  is  given  by: 

S(xa)  =  ikp(x0)\.  (12) 

This  step  is  important  since  an  ideal  dipole  is  not  easily  generated  experimentally.  That  is,  equa¬ 
tion  ( 12)  allows  the  acoustic  source  to  take  the  form  of  a  monopole  which  can  be  purchased  or 
constructed  from  a  simple  speaker. 

The  final  step  in  formulating  the  reciprocity  experiment  for  dipoles  is  the  measurement  of  the  par¬ 
ticle  velocity  u(xs).  The  present  method  used  a  two  microphone  probe  similar  to  those  typically 
used  for  sound  intensity  measurements.  Specifically,  the  acoustic  particle  velocity  was  estimated 
using  a  discrete  representation  of  the  linearized  Euler  equation: 


1  Ap, 

tcop  Ax 


(13) 
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where  the  pressure  gradient  was  measured  with  a  closely  spaced  pair  of  miniature  microphones. 
The  vector  orientation  of  the  original  dipole  source  at*,  is  preserved  by  positioning  of  the  micro 
phone  pair. 

The  magnitude  of  G ^  can  then  be  calculated  by  combining  equations  8  and  10  through  13: 


1  S/)(Xs) 

S(x0)  Ax 


(14) 


The  implementation  of  dipole  reciprocity  therefore  requires  a  monopole  sound  source  with  known 
magnitude  placed  at  jc0,  and  a  dual-microphone  probe  located  at  xs. 


4.2  Experimental  technique 

A  probe  was  designed  and  constructed  that  allowed  for  the  simultaneous  measurement  of  both  G 
and  Gf).  The  probe  and  its  support  were  designed  to  be  small  enough  to  be  easily  manipulated  and 
traversed  in  a  given  experimental  apparatus  without  creating  significant  changes  in  the  local 
acoustic  field.  Additionally,  the  probe  was  designed  to  hold  two  microphones  at  a  fixed  separation 
( A.v)  in  order  to  recover  an  estimate  of  the  pressure  gradient  Ap(xs)/Ax .  The  orientation  of  these 
two  sensors  needed  to  be  articulated  to  match  a  specified  dipole  vector  direction. 

Figure  20  shows  a  picture  of  the  assembled  probe.  The  sensors  were  Knowles  FG  series  electret 
condenser  microphones  commonly  used  in  hearing  aids.  The  diaphragm  size  was  0.7mm,  with  a 
nominal  casing  diameter  of  2.56mm.  The  microphones  were  mounted  on  a  U-shaped  articulating 
arm.  The  microphone  spacing  was  fixed  at  At=  1 8.5mm.  The  arm  was  mounted  to  a  support  sting 
constructed  from  6,35mm  diameter  aluminum  rod.  The  ability  to  rotate  the  sensor  pair  with  two 
degrees  of  freedom  allowed  for  alignment  with  a  given  direction. 

Calibration  of  the  two  miniature  microphones  utilized  a  plane  wave  tube  constructed  from  a 
0.75m  length  of  extruded  PVC  pipe  with  an  inner  diameter  of  40.4mm.  The  cut-on  frequency  for 
non-planar  waves  was  4.87kUz.  A  calibrated  ACO  Pacific,  Inc.  1/4”  free-field  microphone 
(model  7016PM  with  4012  pre-amplifier)  was  used  to  measure  the  plane  wave  amplitude  at  the 
same  axial  location  as  the  reciprocity  probe.  Broadband  acoustic  excitation  was  provided  by  a 
small  speaker  located  at  the  end  of  the  tube.  Analog  signals  were  converted  to  digital  voltage  time 
series  and  recorded  for  post-processing.  The  magnitude  and  phase  of  typical  calibration  functions 
for  the  two  microphones,  //|(ca)  and  H2((o)  are  shown  in  Figure  21.  The  microphone  pair  was 
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observed  to  be  relatively  phase  matched  over  the  frequency  range  shown.  The  average  phase  mis¬ 
match  was  0.05  radians. 

The  measurement  of  G(xs,  xQ)  was  executed  as  follows.  A  B&K  omnisource  (Type  4295)  speaker 
was  positioned  at  the  desired  xG  location.  The  miniature  probe  was  located  at  xs.  The  same  refer¬ 
ence  microphone  used  during  calibration  was  placed  in  the  free  field  which  allowed  for  the  deter¬ 
mination  of  the  source  strength  using  the  free  space  relation: 

|S(o)l  =  (<o)Aa>  (15) 


where  R  is  the  radial  distance  from  the  reference  microphone  to  the  acoustic  source,  and  Ap  is 

the  autospeclral  density  of  the  reference  pressure1,  p,^.  G{xsx0)  was  calculated  by  combining 
equations  4.  9  and  15: 

=  rJ^/  (,6) 


where  the  frequency  band- width,  Ato,  used  in  the  spectral  processing  cancels  so  long  as  they  are 
the  same  for  the  evaluation  of  the  two  spectral  densities. 

Similarly,  equations  14  and  15  can  be  combined  to  obtain  a  simple  form  for  the  dipole  function: 


'Y  Prtf 

where  the  autospectral  density  of  the  pressure  difference.  was  obtained  from 


(17) 


=  Api  +  Api-2Re(ApiPi)  (18) 

and  Ap  jl  is  the  cross  spectral  density  of  the  microphone  pressures. 

The  uncertainty  in  the  measurement  of  G(xs ,  jt„)  was  determined  by  estimating  potential  errors  in 
the  reference  microphone  and  the  probe  microphone  calibrations.  For  the  present  measurements, 
the  reference  microphone  was  calibrated  in  a  standard  B&K.  piston  phone.  The  resulting  error  was 
estimated  to  be  less  than  0.2dB  based  on  the  calibration  uncertainty.  Additionally,  if  the  experi- 


1.  The  definition  and  calculation  algorithms  for  the  aulospectra  density  functions  arc  given  by  Bendat  and 
Piersol  (2000),  Note  that  the  symbol  Gpp  is  commonly  used  for  auto  spectral  density.  However,  G  is 
reserved  here  for  the  acoustic  Cirecn’s  function*  and  the  notation  Ap  will  be  used  to  represent  spectral 
quantities. 
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ments  are  not  conducted  in  a  fully  anechoic  environment,  the  use  of  the  free  space  Green’s  func¬ 
tion  in  equation  15  will  result  in  a  bias  error  in  the  acoustic  source  magnitude.  Issues  related  to  the 
assumption  of  anechoic  boundary  conditions  will  be  discussed  in  the  following  section. 
Uncertainty  in  the  measurement  of  GD  is  related  to  the  measurement  of  the  particle  velocity  using 
the  two-microphone  technique.  This  is  identical  to  the  issues  considered  for  common  intensity 
measurements,  and  detailed  uncertainty  analysis  can  be  found  in  (Seybert  1981,  Fahy  1989,  Wong 
1991.  Jacobsen  &  de  Bree  2005).  For  the  present  investigation,  the  uncertainty  was  found  to  be 
controlled  by  two  factors.  At  low  frequency,  any  phase  mismatch  in  the  sensors  that  was  not  prop¬ 
erly  accounted  for  in  calibration  leads  to  large  errors.  At  high  frequency,  the  approximation 

^  ^  becomes  invalid.  This  type  of  error  is  referred  to  as  systematic  error  by  Fahy  ( 1 989). 

Other  sources  of  error  such  as  electronic  noise  and  calibration  amplitude  factors  were  found  to  be 
substantially  smaller  than  the  phase  mismatch  and  systematic  errors. 

The  phase  error  introduced  when  measuring  the  acoustic  particle  velocity  can  be  found  by  starting 
from  Euler's  equation  for  a  harmonic  pressure  field: 


(19) 


where  $  is  the  true  phase  difference  in  the  pressure  at  two  points  spaced  At  apart.  A  normalized 
uncertainty  can  be  found  by  dividing  by  the  plane-wave  reference  values  and  taking  the  partial 
derivative  with  respect  to  <j> : 


(20) 


where  5<b  is  the  phase  mismatch,  and  the  phase  <|)  =  kAx .  This  error  is  shown  in  Figure  22  using 
S(J>  =  0.05  radians.  This  value  was  considered  to  be  a  conservative  estimate  of  the  phase  uncer¬ 
tainty  based  on  repeatability  and  scatter  in  the  calibration  measurements. 

The  normalized  systematic  error  related  to  the  finite  difference  approximation  of  dp/d. x  is  given 
by  Fahy  (1989)  as 

_  sin(*Av/2) 

*y*  kAx/ 2  1  J 

and  is  also  shown  in  Figure  22, 
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The  uncertainty  defined  by  equations  20  and  21  clearly  define  a  range  of  k&x  for  which  the  mea¬ 
surement  technique  will  be  viable.  A  maximum  normalized  error  of  0.6  (2dB)  limits  the  technique 
to  the  range  0.057  <  kAx  <  4.3 .  The  microphone  spacing  used  in  the  present  measurements  is 
Ax= 18.5mm,  which  corresponds  to  a  frequency  range  of  170 Hz<f<  12.6 kHz.  However,  as 
noted  earlier,  the  cut-on  frequency  of  the  plane-wave  tube  used  for  calibration  was  4.87 kHz .  The 
true  usable  frequency  range  of  the  probe  for  the  measurement  of  Gp  was  then 

1 70//2  < /<  4,81  kHz . 

4.3  Free  space  measurements 

Measurements  of  G  and  Gp  were  obtained  in  an  anechoic  environment  in  order  to  provide  a  direct 
assessment  of  the  technique  and  its  uncertainty.  The  data  were  collected  in  an  anechoic  chamber 
with  dimensions  9.1m  by  7.3m  by  3.6m.  Fiberglass  wedges  with  0.6m  depth  provided  acoustic 
absorption  at  frequencies  above  100Hz.  A  B&K  omnisource  (Type  4295)  was  placed  in  the  room 
and  its  location  designated  as*,,.  The  location  of  the  reciprocity  probe,  xs,  was  1.7m  from  the 
source.  The  probe  axis  was  aligned  radially  from  the  source.  That  is,  the  vector  direction  defined 

by  the  two  microphones  was  parallel  to  & .  The  omnisource  speaker  was  driven  by  a  white  noise 
signal  generated  using  a  Stanford  Research  Systems  function  generator  (model  DS360)  and 
amplified  through  a  B&K  power  amplifier  (Type  2716).  The  reference  microphone  was  located 
i  .3m  from  the  source.  The  sound  pressure  level  at  the  xs  location  was  approximately  68dB.  Time 
series  data  were  recorded  from  the  two  probe  microphones  simultaneously  with  the  reference 
microphone  and  saved  for  post  processing. 

As  noted  earlier,  reflections  that  result  from  incomplete  absorption  of  the  sound  by  the  fiberglass 
wedges  can  lead  to  two  types  of  bias  error  in  the  measured  transfer  functions.  First,  reflected 
sound  that  is  incident  on  the  ,vy  location  will  result  in  a  positive  bias  in  both  G  and  Gp  when  com¬ 
pared  to  the  free  space  conditions.  Similarly,  reflected  sound  that  is  incident  on  the pirj  location 
will  result  in  a  negative  bias  in  G  and  Gp.  Initial  measurements  indicated  that  these  room  reflec¬ 
tions  would  result  in  2  to  5  dB  of  error.  To  cancel  the  effect  of  these  reflections,  the  measurements 
described  above  were  repeated  at  five  dilferent  locations  within  the  anechoic  room,  and  the  result¬ 
ing  functions  were  averaged. 
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The  analytic  solutions  for  G  and  GD  in  free  space  are  given  by 


and 


G(xa,  -v.) 


X,  | 


(22) 


^cosfB)? 


ik\xa  -  ,r,| 


4Wo-xs\ 


(23) 


respectively  (Pierce,  1989).  Note  that  0  is  the  angle  between  the  wavenumber  vector  and  the  line 
segment  defined  by  the  two  probe  microphones.  As  described  before,  0  =  0  degrees  was  used 
for  the  measurement.  The  free  space  functions  given  in  equations  22  and  23  are  shown  in  Figure 
23  along  with  the  experimentally  determined  values.  The  values  of  G  and  Gp  are  shown  on  a  dB 

scale  defined  by  IOIog|(1(G)  referenced  to  Im'1  and  lOiog^Gy)  referenced  to  1m‘\  respectively. 
The  difference  between  the  experimental  and  analytical  values  was  found  to  be  within  +/-  3d B, 
The  bias  in  the  monopole  measurement  over  the  valid  frequency  range  was  0.4dB  with  a  maxi¬ 
mum  error  of-3dB  measured  at  4.32kHz.  The  average  bias  in  the  dipole  measurement  was  0.5dB 
over  the  valid  frequency  range.  The  bias  was  0.7dB  over  the  frequency  range  l00-400Hz.  The 
maximum  error  was  3dB  measured  at  120Hz. 


4.4  Application  to  finite  duct  acoustics 

The  motivation  for  the  development  of  the  present  experimental  technique  was  the  study  ofaeroa- 
coustic  sources  in  low  speed  rotors  operating  in  a  short  duct  section.  The  study  of  the  radiated 
sound  from  such  devices  is  complicated  by  the  presence  of  end  reflections,  organ  pipe  resonance, 
cut-on  of  higher  order  modes,  and  the  acoustic  interaction  of  the  duct  ends  from  outside  the  duct. 
It  is  therefore  of  interest  to  present  initial  results  from  these  studies  to  highlight  the  utility  of  the 
developed  measurement  technique. 

A  schematic  representation  of  the  duct  is  shown  in  Figure  24.  The  duct  walls  were  constructed 
from  extruded  PVC  with  a  machined  inner  radius  a  =  0. 103m  and  a  wall  thickness  of  6.35mm. 
The  length  of  the  duct  was  6a.  A  partial  flange  with  25.4mm  outer  diameter  was  placed  at  both 
ends  of  the  duct  as  shown.  The  location  ofxs  was  centered  both  axially  and  radially  in  the  duct, 
and  the  vector  orientation  was  aligned  with  the  duct  axis.  The  location  of.v0  was  1 .3 1  m  from  the 
duct  center  along  the  duct  centerline. 
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The  measured  functions  G  and  GD  are  shown  as  a  function  of  frequency  in  Figure  25,  The  fre¬ 
quency  was  normalized  by  the  organ-pipe  resonance  frequency: 

fo  “  jf"  (24) 

lLeff 

where  Lejy  =  L  +  2  ■  0.6!  I  a  is  the  effective  length  of  the  duct  using  the  end  correction  factor  of 

Levine  and  Schwinger  (1948).  Also,  shown  on  the  upper  abscissa  are  the  corresponding  values  of 
the  nondimensional  wavenumber,  ka.  Several  notable  features  can  be  seen  in  the  resulting  transfer 
functions.  First,  both  G  and  Gq  are  approximately  the  same  overall  magnitude  as  the  free-space 
functions  given  by  equations  22  and  23  (also  shown  in  Figure  25).  Also,  local  maxima  are  noted 
in  G(oj)  at  odd  multiples  of f„  whereas  G/j(to)  was  found  to  exhibit  these  maxima  at  even  multi¬ 
ples  of/^.  This  result  is  intuitive  based  on  the  symmetry  of  the  geometry  and  the  source  location. 
At  higher  frequencies,  a  number  of  resonances  and  anti-resonances  are  noted.  The  numbered  pairs 
noted  on  the  figure  correspond  to  the  cut-on  frequency  of  higher  order  duct  modes. 

4.5  Conclusions  of  Reciprocity  Measurements 

The  acoustic  transfer  functions  for  single  point  monopole  and  dipole  sources  can  be  represented 
by  the  acoustic  Green's  function  and  its  gradient,  respectively.  The  experimental  determination  of 
these  functions  using  a  reciprocal  technique  has  been  demonstrated.  This  required  the  construc¬ 
tion  of  a  miniature  dual-microphone  probe  that  could  be  positioned  at  a  desired  source  location. 

A  basic  uncertainty  analysis  as  well  as  measurements  conducted  in  a  nominally  anechoic  environ¬ 
ment  have  shown  that  the  technique  provides  an  accurate  representation  of  the  acoustic  transfer 
functions,  w  ith  a  maximum  error  of  about  3dR. 

The  motivation  for  developing  the  present  technique  was  to  provide  a  simple  method  for  the 
determination  of  G  and  GD  in  complex  geometries  for  which  analytical  and  numerical  solutions 
are  either  difficult  to  obtain  or  not  available.  A  demonstration  of  the  technique  using  a  short,  rigid 
duct  section  was  provided.  Experiments  are  presently  underway  for  investigating  the  dependence 
of  Gd  on  the  source  location  for  complex  rigid  and  elastic  duct  sections. 
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Figure  19  :  Illustration  of  acoustic  reciprocity  applied  in  fluids.  The  example  shows  monopole 
(top)  and  dipole  (bottom)  point  sources  located  in  a  cylindrical  duct,  with  an  observation  point 

outside  the  duct. 


Figure  20  :  Probe  used  in  the  measurement  of  the  acoustic  transfer  function. 
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Figure  21  :  Magnitude  (upper)  and  phase  (lower) 
of  the  calibration  transfer  functions  for  both  probe  microphones. 


Figure  22  :  Phase  mismatch  (solid)  and  systematic  (dashed)  normalized  error 
in  the  measurement  of  acoustic  particle  velocity  as  a  function  of  nondimensional  wavenumber. 
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Figure  23  :  Measurement  of  the  monopole  (lower)  and  dipole  (upper)  free  space  acoustic  transfer 
function.  Results  presented  as  t  Olog  )0(G)  (re  1m'1)  and  IO!og|0(GD)  (re  Jrn'2),  respectively. 
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Figure  24  :  Schematic  illustrating  the  cylindrical  rigid  wailed  duct. 
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Figure  25  :  Acoustic  transfer  function  (in  dB)  between  a  monopole  (solid)  and  dipole  (dotted) 
located  at  the  center  of  the  duct  and  an  observation  point  outside  the  duct  on  centerline*  Smooth 

curves  represent  free-space  solutions* 
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5.0  DETERMINATION  OF  DIPOLE  SOURCE  FROM  SOUND  DATA 

This  section  describes  a  method  for  analyzing  the  sound  radiated  by  a  low-speed  axial  How 
ducted  fan.  Specifically,  an  algorithm  was  developed  that  allows  the  determination  of  aeroacous- 
tic  source  strength  from  microphone  measurements.  The  radiated  sound  depends  on  both  the 
sound  generated  by  the  rotor  and  the  geometry  of  the  surfaces  surrounding  the  source  region. 
These  surfaces  act  to  modify  the  spectral  content  of  the  sound  as  perceived  by  an  observer  outside 
of  the  duct.  The  fan  was  considered  as  a  compact  region  of  axial  dipole  sources,  so  that  the  radi¬ 
ated  sound  could  be  described  as  the  product  of  a  frequency-dependent  unsteady  force  and  an 
effective  transfer  function. 

The  microphone  measurements  presented  in  this  section  were  acquired  from  the  radiated  sound 
generated  by  a  ten-b laded  fan  operating  inside  a  finite  length  rigid  duct.  Two  fundamental  fre¬ 
quencies  were  found  to  be  inherent  in  the  system,  one  based  on  the  parameters  of  the  sound 
source,  and  the  other  based  on  the  geometry  of  the  duct.  As  previous  studies  have  found,  these 
two  frequency  scales  allow  for  the  separation  of  the  radiated  sound  into  source  and  transfer  func¬ 
tion  spectra.  A  unique  method  for  determining  these  separate  functions  is  presented,  with  the  pri¬ 
mary  assumption  that  the  aeroacoustic  sources  follow  a  power  law  relationship  with  rotor  tip 
relative  speed.  The  resulting  algorithm  was  shown  to  effectively  quantify  the  acoustic  transfer 
function.  A  set  of  mathematically  defined  input  functions  was  used  to  evaluate  the  efficacy  of  the 
separation  method.  Additionally,  the  ducted  rotor  system  was  operated  over  a  range  of  flow  rates 
to  produce  different  sound  sources,  and  the  transfer  function  calculated  by  the  algorithm  was 
found  to  be  independent  of  this  source  mechanism. 

5.1  Background 

Axial  flow  ducted  rotors  are  important  in  numerous  engineering  applications,  ranging  from  venti¬ 
lation  equipment  to  jet  engines.  The  aerodynamic  sound  radiated  from  these  devices  is  generated 
by  a  number  of  different  mechanisms.  At  low  speeds  these  mechanisms  can  be  generally 
described  by  the  response  of  the  rotor  to  turbulent  flows.  For  example,  the  interaction  between  the 
blades  and  turbulence  in  the  approach  (low  is  often  referred  to  as  “inflow  noise."  Sound  generated 
by  the  scattering  of  unsteady  pressure  that  is  created  by  the  blades  is  often  referred  to  as  “self 
noise."  Many  of  these  sound  sources  have  received  extensive  analytical  attention,  much  of  which 
can  be  found  in  books  by  Blake  (1986)  and  Goldstein  (1975).  and  review  papers  by  Morfey 
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(1973),  Wright  (1976)  and  I  luff  (1997).  Some  examples  of  more  recent  work  include  rotor  inter¬ 
actions  with  the  duct  boundary  layer  (Moissev  et  al  (1978)),  tip  vortex  noise  (Dunne  and  Howe 
(1997)),  vortex  sheet  noise  (Lee  et  al,  (1993))  and  ingested  flow  irregularities  (Atassi  et  al. 
(1993))  and  Kobayashi  (1978). 

The  propagation  of  acoustic  waves  from  rotating  sources  can  be  represented  by  an  appropriate 
Green's  function.  Sound  propagation  from  rotors  operating  in  a  free-space  environment  can  be 
found  by  using  a  free-space  Green's  function  with  modifications  for  rotational  effects  (see.  e.g., 
Blake,  (1986)).  Rotors  operating  in  an  effectively  infinite  duct  invoke  significantly  more  complex 
phenomena  due  to  the  effects  of  the  source  distribution  on  the  excitation  of  various  duct  character¬ 
istic  modes  (see  Doak  ( 1 973),  Pierce  (1981).  and  others).  Noise  that  is  radiated  from  rotors  oper¬ 
ating  in  a  relatively  short  duct  section  exhibits  additional  complexity.  For  example,  reflected 
waves  from  the  duct  ends  will  result  in  "organ-pipe”  modes  whose  amplitude  will  depend  on  the 
location  of  the  sources.  The  directivity  and  amplitude  of  sound  radiated  from  each  of  the  duct 
ends  depends  on  the  specific  geometry  of  the  duct  inlet  and  exit  sections.  The  Green's  function  for 
a  finite  length  duct  can  not  be  determined  analytically,  and  so  must  be  measured  or  computed 
numerically.  For  example,  Boundary  Element  Methods  (BEM)  are  effective,  but  require  signifi¬ 
cant  computational  resources,  particularly  for  very  complex  geometries. 

The  relationship  between  the  generated  sound  and  the  far  field  sound  propagation  can  be  simpli¬ 
fied  significantly  by  considering  the  sound  source  to  be  a  net  equivalent  dipole  force,  F(St,  Mtj  ) . 
The  tip  relative  Mach  number  Mtjp  is  de lined  as  the  relative  fluid  velocity  at  the  tip.  Utif).  divided 
by  the  speed  of  sound  in  the  working  fluid,  c0.  The  independent  frequency  variable  is  described  as 
a  Strouhal  number,  in  the  manner  of  Wcidemann.  and  defined  as 


/Or 

hip 


where  Dr  is  the  rotor  diameter. 


l 

UB  BPF' 


(25) 

Vtjp  is  the  velocity  of  the  rotor  tip, /is  the  dimensional  frequency 


and  li  is  the  number  of  rotor  blades.  For  this  definition  the  fundamental  frequency  of  the  sound 
source  is  taken  to  be  the  rotor  blade  passing  frequency  (BPF).  The  effects  of  the  Green's  function 
can  be  considered  as  a  net  transfer  function.  T(He)s  where  the  Helmholtz  number  is  defined  as 
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„  fD  „  w  B 

He  —  —  =  x  St  x  — T 

CO  JT 

(26) 

where  £>  is  the  duct  diameter. 

The  sound  spectrum  at  a  point  in  the  far  field  is  then  given  by  the  product,  P=F*T.  A  complete 
derivation  of  this  relationship  and  its  assumptions  is  given  in  the  following  section.  The  focus  of 
the  present  investigation  is  the  use  of  experimental  measurements  of  the  radiated  sound  spectra, 
P(f,Mtip),  from  a  variable  speed  ducted  rotor  to  quantify  the  two  functions  FfStM,^  and  T(He). 
These  measurements  utilized  a  tow  speed  fan  operated  inside  of  a  short  annular  duct  section 
inside  an  anechoic  room.  The  radiated  sound  spectra  were  acquired  using  a  single  far  field  micro¬ 
phone  at  1 1  rotational  speeds. 

A  number  of  previous  investigations  have  explored  methods  for  separating  source  and  transfer 
effects  from  measurements  of  the  radiated  sound,  including  Weidemann,  Neise  (1975).  Bent  and 
McLaughlin(1993),  Mongeau,  Thompson  and  McLaughin  (1995)  and  Zhang.  Mongeau  and 
Frankei  (2002).  A  thorough  discussion  of  these  and  other  related  papers  will  be  given  in  the  fol¬ 
lowing  section,  along  with  the  relevant  background  material.  A  new  code  was  written  for  each  of 
these  methods  and  applied  to  the  present  data  set.  All  of  these  methods  were  found  to  identify  the 
majority  of  the  spectral  features  in  both  F  and  T.  However,  in  each  case  the  results  were  found  to 
exhibit  spurious  trends  in  both  T and  F.  Specifically,  these  methods  were  found  to  result  in  func¬ 
tions  for  T  that  were  strongly  dependent  on  the  spectral  content  of  the  sound  source.  It  will  be 
show'n  in  this  paper  that  these  errors  are  likely  due  to  a  non-uniqueness  issue  that  has  not  previ¬ 
ously  been  described  and  accounted  for.  Specifically,  it  will  be  shown  that  large  scale  trends  in 
either  F  or  7’ can  not  be  determined  by  simply  vary  ing  the  fan  speed. 

The  main  result  of  this  section  is  the  presentation  of  a  new  formulation  and  algorithm  for  deter¬ 
mining  F’and  T  from  the  available  measurements.  The  background  and  assumptions  are  provided 
in  the  following  section.  The  efficacy  of  the  algorithm  was  tested  using  a  new  experimental  data 
base  as  well  as  mathematically  defined  input  functions.  The  experiments  included  measurements 
with  dramatically  varied  acoustic  source  spectra  that  resulted  from  throttling  the  fan  into  rotating 
stall.  The  functions  T  returned  from  these  tests  varied  by  less  than  0.29  dR  across  most  of  the  fre¬ 
quency  range  of  interest. 
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5.2  Analytical  background 

5.2.1  Derivation  of  net  equivalent  product  model 

A  system  that  can  be  characterized  as  a  compact  distribution  of  aeroacoustic  sources  surrounded 
by  a  specified  geometry  gives  rise  to  the  conceptual  framework  of  a  net  equivalent  point  source 
and  an  effective  acoustic  transfer  function.  This  concept  can  be  formalized  by  starting  with  the 
Helmholtz  equation, 

Vap  +  =  v  ■  F 

(27) 

where  p  is  the  Fourier  transform  of  the  acoustic  pressure  and  k  is  the  acoustic  wave  number.  1  he 
form  of  the  acoustic  source  on  the  right  hand  side  is  assumed  to  be  dipole  in  nature  w  here  F  is  an 
equivalent  body  force  distribution.  Eq.  (27)  is  valid  for  low  Mach  number  systems  without  mono¬ 
pole  or  quadrupole  type  sources. 

The  Helmholtz  equation  can  be  solved  using  a  Green's  function  approach,  where  Green’s  function 
gfr.rj  is  explicitly  defined  by 


V3tf(r,  to)  +  ^Vr,  rn)  =  -<S(r-  r0), 


(28) 


and  the  boundary  conditions  appropriate  to  the  system.  For  Eq.  (28),  r0  is  the  source  location,  r  is 
the  observer  location,  and  8  is  the  Kronecker  delta  function.  The  solution  to  Eq.  (28)  can  be  writ¬ 
ten  as  a  convolution  of  the  Green's  function  over  a  volume  V0  with  the  distributed  sources  V  •  F , 
giving 


(29) 


where  the  subscript  on  the  operator  V  denotes  operation  over  the  source  region.  A  detailed  deri 
vation  of  Eq.  (29)  is  available  in  Morse  and  Ingard  (1968),  among  others.  Integration  by  parts 
allows  the  far  field  solution  to  be  rewritten  as 


(30) 


If  the  source  region  is  sufficiently  compact  the  Green's  function  is  separable  into  near  and  far  field 
components.  By  using  a  coordinate  system  centered  about  the  source  region,  this  assumption  of 
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compactness  is  denoted  as  kr0< I .  For  the  present  purposes,  it  will  be  useful  to  describe  the 
Green's  function  as  a  product  of  three  functions, 


9  (r.  r0)  =  X  (r)  yfa  (r,  r0)  gnear  (r0) , 


where  gncar  is  a  dimensionless  function  of  the  source  location  that  quantifies  near- fie  Id  effects 
such  as  Doppler  shifts  due  to  moving  sources,  geometrical  non-compactness,  and  other  effects 
due  to  both  propagating  and  non-propagating  components.  The  dimensions  of  the  Green's  func¬ 
tion  are  retained  in  the  free-space  Green's  function.  gfs  =*  — j^_-r--e'fc|r~roL  The  function  x 


accounts  for  amplifications  and  attenuations  at  particular  frequencies  due  to  the  boundary  condi¬ 
tions  imposed  by  the  rigid  duct.  A  similar  amplification  function  was  introduced  by  Morse  and 
Ingard  (p.  320).  As  an  example,  if  a  harmonic  monopole  source  were  placed  near  an  infinite  rigid 
wall,  the  magnitude  amplification  function  would  be  |%|  =  2,  indicating  a  doubling  of  the  magni¬ 
tude  of  the  free-space  Green’s  function  across  all  frequencies. 


Combining  Eq.  (30)  with  Eq.  (3 1)  gives 

P  (r)  =  \  (r)  fff  F  ■  Vo  \gfa  (r.  r0)  gnear  (ro)]  <%. 

JlJ  (32) 

For  the  ducted  rotor  system  considered  in  this  paper,  a  radiation  efficiency  term  can  be  partially 

derived  by  considering  the  acoustic  sources  to  be  axially  oriented  dipoles.  This  is  appropriate  so 
long  as  the  stagger  of  the  blades  is  such  that  the  unsteady  forcing  on  the  rotor  blades  acts  prima¬ 
rily  in  the  axial  direction.  An  equivalent  source  distribution  can  be  defined  as 


F  (r0)  =  Fr  (r„) , 

where  the  unsteady  forces  act  only  in  thex-dircction,  that  is  along  the  axis  of  the  duct,  denoted 
with  the  subscript  x.  The  gradient  operator  in  Eq.  (32)  can  be  carried  out  giving 


(33) 


/'(r)  =  \  (r)  UJ  F,  ■  [V^,  (r.rn),v„,ar  +  gj,  (r,r0)  Vng„riir  (r0)]<JVo- 


I  he  gradient  of  the  free  space  Green's  function  in  the  far  field  is  given  as 


(34) 
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(35) 


Vofir/,  (r,  r0)  =  -  [ikco&e]  gft  (r,  r0) , 

w  here  9  is  the  angle  between  the  axis  of  the  dipole  forces  and  the  vector  between  the  source  and 
the  receiver.  r-r0.  Substituting  this  in  Eq.  (34)  gives 


Pi r)  =  (r)  {ik)9h  (r)  oo *6  fff  F,  -  \ynfar  (r0)  + 

L 


dV  o. 


(36) 


If  r-ra  can  be  approximated  as  |r|,  then  gfs  can  be  removed  from  the  integral  above,  which  is 
appropriate  now  that  the  gradient  operator  has  been  evaluated  and  the  dependence  ofgy;,  on  ra  is 
small.  This  gives 


V  (»')  =  v  (r)  (r)  cos  0  JjJ . 


(ro)  4- 


T* M  (i'o)  1  if  - 

- 1 - 5 — 

ik  cos  0 


The  auto  spectrum  of  the  acoustic  pressure  can  be  obtained  by  taking  the  expected  value  of  Eq. 
(37)  multiplied  by  its  complex  conjugate  to  obtain 

V  -  ¥ (r)  =  l\  (r)[3  A’2  \9f.  (r)|2  cos2  6  fjf  ■  JJj  F^F% 

Vi  Vi 

^019  fiteir  (rDl}l  L  *  v  t  ^Q29»eur  (^ih) 


(37) 


(  *  0 1 ) 


ik  COS  (9 


J  L 


9nrar  (*0?)  A~ 


ik  cm  0 


(38) 

Similar  forms  of  this  expression  have  been  derived  previously  by  Blake  (1986)  and  Zhang  et  ah 
(2002),  among  others.  The  derivation  presented  here  is  unique,  however,  due  to  the  inclusion  of 

the  k“  efficiency  term  and  the  function  x  ■  The  net  equivalent  product  model  can  now  be  given  by 
representing  the  spatial  integral  of  the  two  point  correlation  of  the  distributed  sources  as  a  single 
source  function  F.  yielding 


T  =  TT. 


(39) 


The  use  of  the  notation  T  for  the  acoustic  transfer  function  is  unique  to  the  present  communica¬ 
tion.  and  is  used  to  denote  the  combined  variables 


T  —  |x  (r)|2  fc2  |.9/„  (r))2  cos2  <9 


(40) 


or 
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log  7  =  log  k2  +  log  \gfs  |2  +  log  lx  I2  +  log  cos2  0. 
which  is  valid  for  a  compact  distribution  of  axial  dipoles. 


(41) 


5.3  Separation  of  source  and  transfer  functions 

The  utility  of  the  linear  source-filter  model  given  by  Eq.  (39)  is  apparent  both  for  analysis  of 
experimental  data  and  for  the  prediction  of  acoustic  signatures  from  complex  systems.  A  partial 
separation  of  F  from  T  given  a  set  of  acoustic  measurements  P  is  possible  because  the  acoustic 
sources  are  explicitly  dependent  on  a  frequency  related  to  the  fluid  velocity,  whereas  the  transfer 
function  is  related  only  to  the  geometry  of  the  surrounding  surfaces  and  the  sound  speed  in  the 
working  fluid.  Equation  (39)  can  be  given  as 

log  P  (/)  =  log F  {St)  +  log 7  {He) ,  ^ 

where  it  is  now  clear  that  Eq.  (42)  has  been  reduced  to  two  terms  and  that  F  is  dependent  on  St  but 
not  on  He,  while  T  is  dependent  on  He  but  not  on  St. 

Before  considering  specific  methodologies  for  separating  the  source  and  transfer  functions,  it  will 
be  useful  to  develop  some  insight  into  the  limitations  of  the  analysis.  Consider  a  velocity  depen¬ 
dent  source  function.  F(St),  which  is  broadband  in  nature,  and  a  transfer  function  with  a  single, 
narrow  band  tone  at  a  fixed  frequency.  An  "organ-pipe"  type  of  resonance,  for  example,  could 
result  in  such  a  feature.  As  the  fan  speed  increases,  the  radiated  sound  pressure  would  likely 
increase  with  rotational  speed  according  to  a  power  law.  However,  each  of  these  measured  spectra 
would  exhibit  the  spectral  feature  due  to  the  pipe  resonance  at  the  fixed  frequency.  Separating  this 
feature  of  the  acoustic  transfer  function  from  the  source  function  is  relatively  straight  forward  in 
this  case.  In  the  converse  of  this  example,  consider  a  featureless  transfer  function,  such  as  the  free 
space  Green’s  function,  and  a  source  spectra  with  a  narrow  band  feature,  such  as  a  blade  rate  tone. 
In  this  case,  varying  the  fan  speed  will  result  in  radiated  sound  spectra  in  which  the  frequency  of 
the  narrow  band  feature  is  proportional  to  the  rotor  speed.  Separating  the  features  of  these  func¬ 
tions  is  again  straight  forward. 

It  is  clear  from  Eq.  (42),  however,  that  absolute  magnitudes  of  F  and  T  can  not  be  uniquely  deter¬ 
mined  using  only  measurements  of/1.  What  is  less  obvious,  is  that  broad  trends  can  not  be  deter¬ 
mined  by  making  finite  variations  in  St  or  He.  For  a  simple  example,  assume  a  measured  sound 
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spectra  which  is  constant  with  frequency.  Such  a  spectra  could  be  generated  by  a  source  spectra 
that  increases  with  frequency  combined  with  a  transfer  function  that  decreases  with  frequency,  or 
vice-versa.  Thus,  an  important  aspect  of  any  algorithm  designed  to  separate  F  from  T  must  delin¬ 
eate  which  features  of  the  acoustic  transfer  function  can  be  determined  unambiguously,  and  which 
features  and  trends  can  not  be  determined. 

Previous  attempts  at  separating  the  acoustic  transfer  function  from  the  acoustic  sources  are 
described  by  Weidemann  and  Neise  (1975)  for  centrifugal  rotors.  A  study  by  Neise  and  Barsikow 
(1982)  used  a  series  of  dimensionally  similar  centrifugal  fans  to  study  acoustic  similarity  laws  and 
reported  that  the  acoustic  transfer  function  spectra  were  strongly  dependent  on  the  rotor  flow 
coefficient  (|>,  defined  as  the  mean  flow  velocity  in  the  duct  normalized  by  rotor  tip  velocity.  In 
terms  of  the  Eq.  (42),  the  term  F  should  be  dependent  on  the  flow  coefficient,  but  the  acoustic 
transfer  function  T should  not  be.  In  the  present  derivation,  changes  in  the  radiation  efficiency 
with  flow  coefficient  are  contained  in  the  near-field  term  gneap  which  is  contained  in  F.  The  find¬ 
ing  of  Neise  and  Barsikow  may  have  been  a  result  of  the  underdetermined  aspect  of  the  formula¬ 
tion  and  their  separation  method.  The  results  presented  in  the  current  work  suggest  that  the 
acoustic  transfer  function  T  is  not  a  function  of  rotor  flow  coefficient. 

Mongeau  et  al.  (1995)  developed  a  method  for  characterizing  the  transfer  function  for  a  ducted 
centrifugal  compressor.  A  similar  technique  was  later  applied  to  a  confined  jet  flow  by  Zhang  et 
at.  (2002).  Their  method  requires  the  assumption  that  on  a  plot  of  non-dimensional  sound  pressure 
spectra  in  terms  of  St,  lines  of  constant  He  are  parallel,  with  an  offset  determined  by  the  derivative 
of  the  transfer  function  with  respect  to  He.  This  assumption  is  valid  so  long  as  the  acoustic  source 
function  F(St)  has  only  linear  variations  with  St.  For  example,  if  one  line  of  He  overlaps  a  blade 
rate  tone,  and  another  He  does  not,  the  two  lines  would  not  be  parallel.  These  algorithms  were 
found  to  be  effective  at  determining  the  small  scale  features  of  the  spectral  functions.  However, 
application  of  the  method  presented  by  Zhang  et  al.  (2002)  to  the  current  measurements  led  to 
large  scale  trends  in  the  acoustic  transfer  function  that  were  strongly  dependent  on  the  functional 
form  of  the  acoustic  sources. 

The  present  technique  was  developed  such  that  the  small  scale  features  of  Tcould  be  determined 
unambiguously,  and  the  large  scale  trends  could  be  intentionally  fixed  by  utilizing  additional 
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information  known  about  the  acoustic  system.  The  method  assumes  a  function  F0  exists,  such  that 
the  source  spectra  at  all  rotor  speeds  is  given  by 


?(St.Mtip)=Fn(Si)\Cst\ 


or 


(43) 


(44) 


log,^(S<,  A/np)  =  log/o(Si)  4-  n  (St) logMip- 

In  words,  increasing  tip  speed  will  result  in  a  linear  increase  in  the  logarithm  of  generated  sound 
pressure  at  a  fixed  value  of  Si.  The  purpose  for  leaving  n  as  a  function  of  Si  in  Eq.  (44)  was  two 
fold.  First,  n  can  be  found  as  a  function  of  St  by  a  least  squares  fit  to  F  and  Mlip  using  Eq.  (44). 
Thus  leaving  n  as  a  function  of  St  seemed  to  provide  the  least  restrictive  assumptions.  Second,  the 


tip  speed  scaling  P~Vn  is  given  significant  attention  in  the  literature,  and  many  different  values 
for  tonal  and  broadband  sound  have  been  documented;  see.  for  example,  Neise  (1975).  The  meth¬ 
ods  developed  by  Mongeau  et  al.  (1995)  and  Zhang  et  al.  (2002)  require  a  fixed  value  of  n  but 
provide  results  that  are  explicitly  independent  of  the  value  of  n  chosen.  A  study  by  Quinlan  and 
Krane  (1996)  reported  on  the  spectra!  form  of  the  tip  speed  exponent n,  noting  that  the  velocity 
exponent  should  be  a  a  function  ofS!f,  but  they  approximated  it  as  n(j)  by  using  a  small  change  in 
rotor  tip  speeds  (Clmax/£lmin  =1.16)  such  that  the  “smearing”  of  spectral  features  such  as  blade 
rate  tones  over  Si  was  minimized.  No  such  approximation  will  be  required  of  the  current  work.  A 
detailed  study  of  the  resulting  function  n(St)  is  beyond  the  scope  of  this  paper,  however. 

The  main  point  of  this  paper  is  to  describe  a  method  for  using  Eqs.  (41),  (42)  and  (44)  to  find  the 


function  x  ■  The  method  works  by  assuming  that  P  is  not  proportional  to  V”  only  because  logx  is 


nonzero.  That  is,  radiated  sound  at  the  fixed  value  of  St  should  follow  a  Vn  power  law.  Deviations 
from  this  are  assumed  to  be  a  result  ofamplifications  and  attenuations  at  specific  frequencies  1  le, 
quantified  as  the  function  %  (He). 


5.3. I  Ducted  rotor  experiment 

A  review  of  the  rigid-walled  ducted  rotor  and  the  specifics  of  the  acoustic  measurements  is  pro 
vided  in  this  subsection.  The  experimental  model  was  constructed  in  the  anechoic  chamber.  A 
center  cylinder  extended  upstream  of  the  duct  inlet  and  downstream  of  the  outlet  to  create  an 
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annular  flow  passage.  The  inlet  (low  was  conditioned  using  a  single  layer  of  a  "cheesecloth"  type 
thin  fabric  with  negligible  acoustic  impedance  stretched  over  a  cylindrical  wire  frame. 

Radiated  sound  was  measured  using  B&K  6.35  mm  capacitance  based  microphones.  The  micro¬ 
phones  were  positioned  approximately  two  meters  from  the  duct  inlet  at  an  angle  of  9.5  degrees 
from  the  duct  axis.  Measurements  were  acquired  at  i  1  equally  spaced  speeds  from  2500  to  5000 
RPM.  Example  sound  spectra  are  shown  in  Figs.  26  and  27.  A  number  of  distinctive  features  were 
observed  in  the  spectra.  For  example,  multiple  “organ  pipe"  type  resonances  can  be  observed  at 
fixed  frequencies  in  Fig.  26.  The  spectra  are  plotted  as  functions  of  St  in  Fig.  27  and  show  distinct 
features  at  integer  multiples  of  St,  such  as  blade  rate  tones.  Significant  broad  band  sound  was  also 
observed  at  high  frequencies. 

The  radiated  sound  was  found  to  be  dependent  on  the  rotor  flow  coefficient  The  duct  flow  rate 
was  decreased  by  applying  variable  restriction  in  the  form  of  thin  fabric  with  negligible  effect  on 
the  acoustic  field  to  the  duct  exit,  and  increased  by  attaching  the  outlet  of  the  duct  to  the  anechoic 
wind  tunnel  blower. 


figure  26  Measured  sound  spectra  for  ducted  rotor  at  three  rotational  speeds. 
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Figure  27  Measured  sound  spectra  for  ducted  rotor  at  three  rotational  speeds,  frequency 

normalized  as  St 


5.4  Spectral  decomposition  algorithm 

This  section  describes  an  the  details  of  the  algorithm  to  extract  an  acoustic  transfer  function  spec¬ 
trum  from  a  series  of  microphone  measurements  acquired  at  a  range  of  rotor  tip  speeds.  The  linear 
model  given  in  Eq.  (42)  can  be  expanded  by  substituting  Eq.  (44)  for  log  Fand  Eq.  (4 1 )  for  log  T 
to  obtain 

log  ’P(f,  Mnp)  =  log Fb(St)  +  n(St)  log  A/tip 

+  log  k2  +  log  |  <jf,  |3  +  log  |\(//e)| 2  +  log  con3  9. 

In  this  expression  Fa(St).  n{St)  and  %(He)  are  unknown.  The  algorithm  seeks  to  find  the  function 
X  (He)  that,  in  a  least  squares  sense,  best  fits  this  expression  given  measurements  of  P  at  several 
values  of  M,  ip.  Once  the  function  x(He)  has  been  found,  the  value  of  the  true  acoustic  transfer 
function  can  then  be  reconstructed  using  Eq.  (41).  The  relationship  between  the  source  function, 
the  transfer  function,  and  the  radiated  sound  is  given  in  Eq.  (42). 
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5.4.1  Implementation  as  a  computer  algorithm 

The  method  was  implemented  as  an  algorithm  using  Matlab  on  a  desktop  computer,  Processing 
time  ranged  from  less  than  a  minute  to  several  minutes,  depending  primarily  on  the  number  of 
points  in  the  frequency  spectra.  A  flowchart  representation  of  the  algorithm  is  included  in  Rig.  28 
for  reference.  The  algorithm  is  iterative,  and  the  index j  is  used  in  this  section  to  denote  terms  that 
are  calculated  with  each  iteration. 


Figure  28  Flowchart  representation  of  decomposition  method.  The  number  of  the  iteration 

is  denoted  by  the  index  /'. 


Step  I :  The  source  spectra  for  each  rotor  speed  are  estimated  as 

log Fj(f,  Mt*)  =  log  V{f, Mu,)  -  log*1 

-  log  Iff/, |a  -  log  |  \j{  He)  |J  -  log  cos1  $.  ^ 

For  the  first  iteration  (J  =  1 )  the  function  log  j\;=1  (//<-)!*  is  set  equal  to  zero,  and  the  rest  of  the 

terms  on  the  right-hand  side  of  Eq.  (46)  are  known.  The  impact  of  this  initial  value  for  x  w  ill  be 


t  diversity  of  Noire  Dame 


Center  for  Flaw  Physics  and  Control 


Sound  and  Vibration  Measurements  Using  an  Air-analog  Model 


47 


discussed  below.  The  spectra  are  made  non-dimensional  using  the  source  length  iiDr  and  velocity 
V,ipB  as. 


FjiSt,  A/(,p)  =  I  ^'i,pB/rc Dr )  Tj{  f.  Mhp). 


(47) 

This  type  ot  spectral  scaling  is  described  by,  e.g.  Blake(l986)  (p.  38),  Linear  interpolation  was 
used  to  obtain  the  spectral  magnitudes  at  discrete  values  of  St  from  Fj(f,  Mlip) .  It  was  also  impor¬ 
tant  to  restrict  the  calculations  to  the  extent  of  the  available  data.  Specifically,  the  highest  fre¬ 
quency  usable  was  determined  by  the  highest  St  of  the  highest  speed.  As  an  example,  the 
maximum  measured  sound  frequency  of  20  kHz  corresponded  to  St=  48  at  2500  RPM,  and  St  = 
24  at  5000  RPM,  so  the  highest  St  that  could  be  used  was  24.  The  limit  of  the  estimates  for 
Fj{St=24)  likewise  correspond  to  frequencies  of  10  kHz  to  20  kHz,  so  the  lower  value  of  10  kHz 
must  be  used  in  re-dimensionalizing  the  functions  back  to  frequency  in  Step  4  below.  The  lower 
limit  of  the  usable  range  was  determined  by  the  usable  frequency  range  of  the  anechoic  facility.  A 
minimum  useful  frequency  of  100  Hz  for  the  present  measurements  corresponded  to  St  =  0.24  for 
2500  RPM  and  St  =  0.12  for  5000  RPM.  The  limiting  usable  value  ofSt  =  0.24  was  used,  which 
corresponded  to  100  Hz  at  2500  RPM  and  200  Hz  at  5000  RPM.  The  frequency  range  of  usable 
data  is  reduced  at  both  ends  by  the  scaling. 


Step  2:  A  least-squares  fit  to  the  estimated  sound  source  spectra  is  used  to  obtain  Hj(St)  and 
F0,{St),  using  the  expression  provided  in  Bq.  44,  and  denoted  here  with  iteration  index j  as 


=  l°g  A;  (St)  +  nj(St)  !°g 


(48) 

According  to  this  expression,  at  each  value  ot  St.  a  plot  of  log  Fj{StMtip )  vs  log  I Mti  should  be  a 
straight  line  with  slope  nj(St).  Deviations  from  the  linear  scaling  are  assumed  to  be  a  result  of 
\  ^  x, .  An  example  from  the  ducted  rotor  experiment  is  shown  in  Fig.  29  at  St  =  1 ,02.  In  the  cur¬ 
rent  experiment,  P(St=  1 .02.Mljp)  spans  a  frequency  range  from  425  to  850  Hz  when  measure¬ 
ments  were  obtained  at  1 1  equally  spaced  speeds  between  2500  and  5000  RPM.  The  values  of 
rij(St)  were  determined  using  an  ordinary  least-squares  fit,  as  well  as  a  robust  linear  fit  method,  as 
described  by  Huber  (1981).  The  robust  fit  was  implemented  using  the  “robustfitO”  command  in 
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Matlab,  which  uses  an  iterative  weighted  adaptation  of  the  least-squares  method  that  reduces  the 
influence  of  outliers  on  the  linear  fit.  In  a  small  number  of  cases,  less  than  1%  for  the  present  data 
set,  the  robust  fit  fails  to  converge,  and  an  ordinary  least-squares  fit  is  used  for  these  values  of  St, 


Figure  29  Experimental  data  showing  linear  fit  to  log  Fj=)(St).  estimated  sound 
source  log  f}(St=  1 .02);  solid  line:  robust  fit;  dashed  line:  least-squares  fit. 


Step  3:  An  improved  estimate  for  log r{St  A/Kp)|,  denoted  log  A/„p)  was  obtained  using  Eq. 
(48),  along  with  nfiSt)  and  log  FOJ{Sf)  from  Step  2.  This  step  can  be  expressed  as. 


log  M„p)  —  log  f0'}(St)  +  rijfSr )  log  Mtip. 

This  new  estimate  of  the  source  spectra  is  essentially  a  smoothed  version  based  on  the  linear  fit, 
and  most  of  the  transfer  function  effects  should  have  been  removed. 


(49) 


Step  4;  An  estimate  for  x  can  now  be  found  for  each  rotor  speed  using, 

bg  |  vJ+1  (/,  AM 2  =  log  P(/,  Mttp)  -  log  f}{f,  Mtip) 

-  log  ft3  ~  log  \<lf.  Ia  -  log  cos3  9, 


(50) 
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as  found  in  Eq.  (45).  The  dimensionless  source  spectra  F'^StMup)  were  used  to  obtain  F’jif,Mnp ), 
and  interpolation  was  again  required  to  obtain  data  at  the  same  points  in / as  the  original  measured 
sound  spectra  for  use  in  Eq.  (50).  The  transfer  function  effects  can  now  be  quantified  as 

|\j+i(  Wf)r.  where  the  overbar  denotes  average  over  rotor  speeds  Mtip. 


Step  5:  The  function  l°K  IXj+if  calculated  in  the  previous  step  can  be  used  in  place  of 
I  \j  ( //e)  j J  to  iterate  from  Step  I. 


Iterating  between  one  and  three  times  was  found  to  provide  a  balance  between  a  converged  solu¬ 
tion  and  a  useful  frequency  range  in  the  result.  The  example  data  from  the  ducted  rotor  in  Fig.  29 
are  shown  in  Fig.  30  after  one  iteration.  This  indicates  the  transfer  function  features  have  been 
largely  captured  in  x  (He),  and  the  linear  scaling  is  seen  to  be  significantly  improved. 


Figure  30  Experimental  data  showing  sound  source  linear  scaling  after  one  iteration  of  the  algo 
rithm.  Transfer  function  effects  have  been  largely  removed. 
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5.4.2  Error  analysis  of  the  algorithm 

A  set  of  defined  input  functions  was  constructed  in  order  to  evaluate  the  effectiveness  of  the  algo¬ 
rithm.  The  set  consisted  of  a  function  for  x  (/fe),  a  function  for  F0(St)  and  a  scaling  n{St).  Spectral 
functions  P(fM,iP)  were  created  using  Eq.  (45)  and  used  as  input  into  the  algorithm.  The  resulting 
function  for  \(Ue)  could  then  be  compared  to  the  input,  and  an  error  computed  by  comparing  the 
calculated  transfer  function  (xc{Hc))to  the  actual  transfer  function  { \A(//r))  as 


J  log|xc(//e)|2-log|xA(He)r 


dlle 


N 


j[log|VA(i*e)|a]  dHe 


(51) 


The  transfer  function  log|\A(ffe)|a  used  for  the  defined  input  function  set  was  a  gauss ian  curve 
with  a  parameterized  frequency  width  Ay  and  center  frequency^,  given  as. 


log|\A(ife)2 


e-ioi  f-U)*/*} 


(52) 

I  his  function  is  shown  in  fig.  3 1.  The  feature  shape  can  be  considered  a  worst  case,  since  the  val¬ 
ues  were  only  positive  and  only  induced  outliers  on  one  side  of  the  correct  linear  slope  in  the 
\o%F(St,Mtip)  vs  log  Mtjp.  A  simple  curve  for  the  source  spectra  was  used: 

logjTo  =  30  cos  (5f/2r>)  +5  cos  (2ltSt  -  n)  e_s,j 


Other  functional  forms  of  the  acoustic  source  w'ere  used,  and  the  results  indicate  that  the  function 
\(Hf)  did  not  depend  on  Fa.  A  scaling  of  n  =  const  was  used  for  these  results,  although  identical 
results  were  obtained  when  n(St)  was  used.  The  experimental  rotor  tip  speeds  were  used  for  Mti.r 
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Figure  32  Error  for  different  values  of  feature  size  ratio  and  speed  ratio  after  three  itera¬ 
tions.  Dark  line:  SR=2,  Med-gray:  SR=3,  Light  gray:  SR=5. 
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A  series  of  transfer  functions  were  tested  over  a  range  of fc  and  Ay-and  the  error  was  computed,  it 
was  found  that  the  error  could  be  described  by  the  functional  relationship 


(53) 


where  SR  is  the  ratio  of  the  maximum  rotor  speed  to  the  minimum  rotor  speed.  Specifically,  the 
error  curves  at  each  SR  collapsed  when  plotted  against  the  feature  size  ratio,  Sj!fc.  This  is  expli¬ 
cable  because  the  range  of f  spanned  at  fixed  St  is  a  linear  function  of  St,  In  this  experiment,  at  St 
=  1  the  11  speeds  span  416  to  833  Hz  while  at  St  =  10  the  range  is  between  4160  and  8330  Hz. 
l  ower  St  values  result  in  frequency  ranges  that  are  smaller,  such  that  if  the  transfer  function 
changes  only  slightly  across  the  frequency  range  it  can  not  be  well  resolved.  The  apparent  size  of 
a  feature  therefore  depends  on  its  center  frequency  "location,”^,  as  well  as  its  frequency  "width." 
A f.  Note  that  this  result  is  consistent  with  the  previously  mentioned  limits  on  separating  Eq.  (42). 
Specifically,  absolute  magnitudes  and  linear  trends  in  the  transfer  function  can  not  be  described 
by  finite  values  of  A f!fc  and  thus  can  not  be  determined. 

The  successful  identification  of  transfer  function  features  is  also  dependent  on  the  speed  ratio 
used  during  measurement.  A  small  speed  ratio  results  in  points  that  are  close  together  in  /at  con¬ 
stant  St  and  only  transfer  function  features  that  cause  changes  in  that  range  of  frequency  can  be 
identified.  A  larger  speed  ratio  allows  larger  features  to  be  identified.  Minimum  and  maximum 
rotor  speeds  are  limited  by  the  experimental  apparatus.  Conversely,  the  limitation  on  using  larger 
speed  ratios  is  important  when  considering  the  scaling  the  spectra  between  St  and  He,  as  previ¬ 
ously  discussed.  This  restricts  the  number  of  iterations  that  can  be  performed,  although  adequate 
results  can  typically  be  obtained  w  ith  only  one  or  two  iterations. 

The  error  functions  defined  by  Eq.  (53)  are  shown  in  Fig.  32.  As  shown,  for  the  set  of  defined 
input  function  tests  using  SR  =  2  and  a  feature  of  size  0.1  an  error  of  0.004  was  calculated  after 
one  iteration  and  0.0008  after  two  iterations.  A  larger  feature  w  ith  size  of  0.2  resulted  in  an  error 
value  of  0.05  after  one  iteration  and  0.03  after  three  iterations.  Small  features  are  well  resolved 
under  all  speed  ratios  while  larger  features  are  poorly  resolved.  A  higher  speed  ratio  leads  to 
lower  error,  especially  tor  larger  feature  size  ratios. 
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5.5  Results  using  experimental  data 

An  effective  test  of  the  algorithm  utilized  data  acquired  from  the  ducted  rotor  experiment  using 
the  same  duct  geometry  with  a  variety  of  rotor  sound  sources.  Manipulating  the  rotor  flow  rate 
provided  a  method  for  varying  the  sound  produced  by  the  rotor  without  introducing  additional 
sound  sources  into  the  system.  Ingested  turbulence,  for  example,  could  also  have  been  varied,  but 
turbulence  generators  may  be  a  significant  source  of  sound.  The  rotor  flow  coefficient  <|)  was  var¬ 
ied  between  rotor  stall  and  zero  pressure  rise.  The  spectra  forfi  =  5000  RPM  are  shown  in  Fig.  33 
at  four  flow  coefficients.  Substantial  changes  in  the  sound  produced  by  the  system  can  be 
observed,  indicating  that  the  source  spectra  was  a  strong  function  of  <|>. 

The  function  \(He)  was  calculated  as  specified  in  the  previous  sections  using  microphone  mea¬ 
surements  acquired  at  the  four  flow  coefficients.  These  results  are  shown  in  Fig.  34  using  two  iter¬ 
ations  of  the  algorithm.  The  frequency  axis  was  normalized  by  the  first  organ  pipe  frequency 
based  on  duct  length  L  calculated  as 

10  2(^+0.61330)’  (54) 

using  the  end  correction  as  calculated  by  Levine  and  Schwinger  (1948).  The  harmonics  of  pre¬ 
dicted  organ  pipe  modes  are  seen  to  correlate  with  peaks  in  the  resultant  transfer  function.  It  is 
also  clear  that  the  function  \  (Hr)  determined  by  the  algorithm  is  essentially  independent  of  the 
acoustic  source.  The  standard  deviation  between  the  four  curves  at  each  frequency  was  calculated, 
and  a  mean  value  of  0.29  dB  was  found. 

Significant  deviations  are  noted  at  only  two  frequencies../^35!  and  )%=8.  At  the  lower  frequency, 
the  deviant  \{Hv)  values  result  from  the  lowest  flow  coefficient.  4»  =  0.17.  The  rotor  is  stalled  at 
this  condition  and  the  reduced  magnitude  was  likely  a  result  of  changes  in  the  coupling  between 
the  dipole  source  distribution  and  the  duct  plane  waves.  The  deviant  \(Up)  values  at ,/%=8 
occurred  at  the  highest  flow  coefficient. /=  0.42.  At  this  condition  the  blades  are  completely 
unloaded  and  again,  a  change  in  the  coupling  between  the  blade  dipoles  and  higher  order  duct 
modes  could  explain  the  variation  in  \(Hr), 

The  final  step  in  the  analysis  was  to  reconstruct  the  transfer  function  T(He)  from  the  determined 
\{H> )  function  using  Eq.  (41).  The  sound  source  function  F could  be  found  using  T and  the  mea¬ 
sured  sound  spectra  P \  as  found  in  Eq.  (42).  An  example  of  these  three  functions  is  shown  in  Fig. 
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35.  As  an  independent  validation  of  Eq.  (41)  and  its  assumptions,  a  reciprocal  technique  was  used 
to  obtain  the  dipole  transfer  function.  Specifically,  a  miniature  intensity  probe  was  constructed, 
calibrated  and  placed  in  the  duct  at  the  axial  position  of  the  rotor  with  the  blades  removed.  The 
dipole  transfer  function  (Vos)  was  then  determined  using  the  method  of  Fahy  (1995). 

The  transfer  function  determined  using  reciprocity  is  shown  in  Fig.  35  as  a  dotted  line.  The  abso¬ 
lute  magnitude  and  trend  of  the  two  curves  support  the  present  interpretation  of  Eq.  (41 )  for  St  < 
2,  corresponding  to  around  1 700  Hz.  Note  that  the  features  determined  from  reciprocity  are 
greater  in  magnitude  by  up  to  1 5  dB  compared  to  the  rotor  transfer  function  T at  integer  multiples 
of/0.  This  result  is  likely  due  to  the  absence  of  the  blades  during  the  reciprocity  experiment.  That 
is,  the  rotor  would  cause  an  effective  impedance  change  in  the  duct  which  would  reduce  the  organ 
pipe  resonances.  In  addition,  the  true  spatial  distribution  of  the  sources  (axially  and  radially)  that 
occurs  in  the  rotor  would  lead  to  reduced  resonant  response. 
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Figure  33  Measured  sound  under  different  How  coefficients  at  5000  RPM. 
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Figure  34  Calculated  transfer  functions  at  varied  flow  coefficients.  Frequency  is  normal 
ized  by  the  first  organ  pipe  mode  frequency. 
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Figure  35  Example  sound  spectra  and  decomposition  results  from  algorithm  showing  sound 
pressure,  P,  unsteady  force,  F,  and  transfer  function.  T. 


5.6  Conclusions 

A  technique  for  processing  aeroacousttc  measurements  was  developed  and  evaluated.  Using  this 
algorithm  the  acoustic  transfer  function  T  could  be  determined.  This  will  allow  for  the  study  of  the 
sound  source  spectra  F  separate  from  acoustic  transfer  function  effects.  The  method  utilized  the 
two  frequency  scales  that  are  inherent  in  many  applications  to  determine  the  functional  form  of  a 
net  equivalent  source  and  a  net  acoustic  transfer  function  from  radiated  sound  spectra  at  various 
fan  speeds.  This  was  implemented  as  a  computer  algorithm,  and  evaluated  with  mathematically 
defined  input  functions,  as  well  as  experimental  results  from  a  low  speed  ducted  rotor.  The  fol¬ 
lowing  conclusions  have  been  drawn  from  this  study: 
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•  Any  separation  of  net  sources  and  net  transfer  function  requires  some  form  of  assumption 
about  either  the  source  or  the  transfer  function  spectra.  In  the  present  algorithm,  the  acoustic 
source  was  assumed  to  increase  with  rotor  speed  to  an  unknown  power.  It  was  found  that  this 
assumption  allowed  for  a  unique  determination  of  which  contained  the  narrow  band 
spectral  features  of  the  acoustic  transfer  function  T. 

•  The  ability  of  the  algorithm  to  determine  spectral  features  of  the  net  transfer  function  was 
quantified.  Specifically,  features  whose  width  to  center  frequency  ratio  \j/fc  <  0. 1  were  found 
to  be  correctly  resolved  with  an  error  of  less  than  1%  for  a  speed  ratio  of  2.  This  result  was 
found  to  be  independent  of  the  spectral  character  of  the  acoustic  source. 

•  The  algorithm  was  applied  to  a  ducted  rotor  model  operating  under  a  range  of  How  coeffi¬ 
cients.  Transfer  functions  were  calculated  for  each  of  these  source  conditions  and  were  found 
to  closely  agree,  with  a  standard  deviation  of  0.29  dB.  The  method  effectively  determines 
transfer  function  features  from  measured  sound  in  a  system  where  a  net  dipole  source  term  is 
an  appropriate  approximation. 

•  The  spectral  character  of  the  ducted  rotor  transfer  function  was  found  to  have  peaks  of  magni¬ 
tude  up  to  10  dB  at  integer  multiples  of  the  first  duct  organ  pipe  mode,  in  addition  to  other  fea¬ 
tures. 

•  The  magnitude  and  large  scale  trend  of  the  transfer  function  can  be  determined  when  a  distri¬ 
bution  of  axially  oriented  dipoles  is  an  appropriate  approximation  for  the  sound  source.  The 
magnitude  and  large  scale  trends  were  verified  using  an  independent  reciprocity  measurement 
of  the  transfer  function. 
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6.0  Sound  and  Vibration  Produced  by  an  Airfoil  Tip  in  a  Turbulent 
Boundary  Layer  Flow  with  an  Clastic  Wall 

6.1  Introduction 

Predicting  the  acoustic  radiation  from  a  finite  elastic  wall  section  subjected  to 
complex  turbulent  forcing  is  a  complicated  problem.  This  is  particularly  true  when  rigid 
surfaces  are  subjected  to  turbulent  flow  in  proximity  to  the  vibrating  structure.  This  type 
of  interaction  occurs  often  in  applications  and  results  in  a  coupling  of  rigid  body  aerody¬ 
namics,  aero-acoustics,  and  structural  acoustics. 

A  flowchart  describing  the  two  boundary  value  problems  considered  in  the  present 
work  is  shown  in  Fig.  36.  The  first  row,  labeled  ’A’,  describes  the  canonical  problem  of  a 
membrane  subjected  to  a  stationary  turbulent  flow  on  one  side.  A  formal  solution  for  the 
structural  motion  can  be  found  by  convolving  the  unsteady  wall  pressure  amplitudes 
(expressed  in  a  wavenumber-frequency  domain)  with  the  structural  Green's  function  1. 
Given  the  structural  motion,  the  radiated  acoustic  pressure  can  be  found  through  a  second 
convolution  of  the  wall  motion  with  the  acoustic  Green's  function  [1  ].  Although  problem 
'A'  can  be  solved  formally,  a  number  of  issues  make  practical  predictions  difficult.  For 


Figure  36:  Flow  chart  showing  describing  each  of  the  physical  mechanisms 
involved  in  the  present  boundary  value  problem. 
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example,  the  unsteady  wall  pressure  may  result  from  external  acoustic  excitation  or  turbu¬ 
lence  which  is  not  stationary  in  either  space  or  time,  or  both.  Secondly,  the  structural 
Green's  function  may  not  be  known  analytically,  and  may  be  difficult  to  resolve  numeri¬ 
cally.  This  is  particularly  true  when  the  wave  motion  of  the  structure  includes  mixed  flex¬ 
ural  and  membrane  waves,  and  when  fluid  loading  effects  are  present  [2-5].  Lastly,  mean 
static  pressure  variations  can  lead  to  localized  membrane  stresses  in  the  structure  that 
modify'  its  response. 

The  second  row  shown  in  Fig.  36  illustrates  a  more  complex  problem  in  which  an 
airfoil  tip  is  placed  in  proximity  to  the  elastic  wall.  In  this  configuration  the  leading  edge 
of  the  airfoil  interacts  with,  and  modifies,  the  approach  turbulence,  and  results  in  a  net 
unsteady  force  on  the  airfoil  surface.  This  unsteady  lift  can  result  in  radiated  dipole  type 
sound  which  can  add  to  the  unsteady  pressure  experienced  by  the  elastic  wall,  in  addition, 
the  movement  of  the  membrane  can  lead  to  changes  in  the  pressure  field  experienced  by 
the  airfoil,  leading  to  a  fully  coupled  problem. 

The  present  investigation  considered  the  conditions  with  and  without  an  airfoil  in 
proximity  to  a  round  Aluminum  membrane.  The  membrane  was  inserted  into  one  side  of  a 
test  section  of  an  anechoic  wind  tunnel.  A  detailed  description  of  the  experimental  pro¬ 
gram  and  data  processing  is  provided  in  the  following  section.  Results  from  each  stage  of 
the  causal  flowchart  shown  in  Fig.  36  are  provided  for  cases  'A1  and  'B1  in  sections  6.3 
through  6.6.  Conclusions  will  then  be  presented  and  discussed  in  the  context  of  the  ongo¬ 
ing  research. 

6.2  Experimental  Setup 
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Wall  Pressure  Sensors 


Figure  37:  Schematic  of  experimental  setup. 


A  schematic  representation  of  the  experiment  is  shown  in  Fig.  37.  The  open  jet 
test  section  was  enclosed  on  three  sides  allowing  the  airfoil  to  be  cantilevered,  with  a  tip 
gap  of  6.35  mm,  into  the  test  section.  Acoustic  pressure  generated  on  both  sides  of  the 
membrane  could  radiate  into  the  anechoic  chamber.  The  flat  plate  airfoil  had  a  span  of  406 
mm,  chord  of  610  mm  and  thickness  of  6.35  mm.  The  large  chord  and  span  provided  a 
force  response  with  only  leading  edge  effects  as  described  in  the  next  section.  The  leading 
edge  was  milled  to  a  5  to  1  ratio  ellipse.  Data  were  acquired  for  four  Mow  speeds  ranging 
from  17.5  to  25  m/s  in  increments  of  2.5  m/s. 
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The  elastic  wall  was  a  circular  membrane  ofQ.  127  mm  thick  Aluminum  foil  with  a 
diameter  of  0.4064  m.  A  uniform,  isotropic  tension  of  219  N/m  was  applied  to  membrane. 
The  tension  was  measured  using  four  force  transducers.  The  membrane  was  clamped  on 
the  entire  circumference. 

Two  turbulence  generating  grids  were  combined  in  order  to  produce  the  unsteady 
approach  How  field.  The  first  grid  was  made  from  12.7  mm  diameter  cylinders  with  76.2 
mm  spacing  while  the  second  was  made  from  3.2  mm  diameter  cylinders  with  25.4  mm 
spacing.  The  purpose  of  two  grids  was  to  augment  the  natural  elastic  wall  boundary'  layer 
and  generate  a  large  range  of  scales  of  turbulent  fluid  motions.  Hot  wire  ancmometry  was 
used  to  determine  the  velocity  spectrum  of  the  approach  flow. 

The  unsteady  lift  response  of  the  airfoil  to  the  turbulent  approach  flow  was  mea¬ 
sured  directly  using  novel  "strip"  surface  pressure  transducers  embedded  in  the  tip  of  the 
airfoil  as  shown  in  Fig.  37.  A  sensor  was  placed  on  either  side  ofthe  airfoil  tip  near  the 
leading  edge.  The  shape  of  the  sensors  provide  a  response  that  is  directly  proportional  to 
the  net  unsteady  lift  experienced  by  the  leading  edge  region  of  the  foil.  Data  were  acquired 
at  a  sampling  frequency  of  40  kHz  for  209  s.  A  detailed  description  of  the  measurement 
technique  may  be  found  in  Slaboch  and  Morris. [6] 

The  wall  pressure  spectrum  ofthe  turbulent  flow  field  was  measured  using  an 
array  of  16  microphones  with  a  diameter  of  0.7  mm  embedded  in  a  rigid  wall  as  shown  in 
Fig.  37.  The  microphones  were  flush  mounted  to  the  surface  and  spaced  1.9  cm  apart  in 
the  transverse  direction  with  8  on  each  side  of  the  airfoil. 

The  elastic  wall  vibrational  response  was  measured  using  a  scanning  laser  Dop¬ 
pler  vibrometry  system  (SLDV)  and  a  separate,  single  point  LDV.  The  SLDV  was  used  to 
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survey  145  points  on  the  membrane,  while  the  single  point  LDV  was  focused  only  at  the 
center  point  of  the  membrane.  A  full  scan  of  the  membrane  was  acquired  for  each  of  four 
flow  speeds  both  with  and  without  the  airfoil  present.  Data  were  acquired  at  40.96  kHz  for 
12.8  s  per  scan  point. 

The  far  field  acoustic  noise  generated  by  the  system  was  recorded  using  a  6.35  mm 
ACO  microphone.  The  microphone  was  placed  2  m  from  the  membrane. 

6.3  Unsteady  l  ift  Response 

The  airfoil  configuration  described  in  Section  6.2  was  chosen  specifically  to  iso¬ 
late  the  leading  edge  response.  That  is,  the  trailing  edge  was  sufficiently  far  downstream 
as  to  not  affect  the  membrane  motion.  The  leading  edge  response  to  unsteady  fluid 
motions  is  described  by  Amiet  17].  The  unsteady  lift  can  be  represented  in  non-dimen¬ 
sional  form  by: 


Figure  38:  Unsteady  lift  response  of  the  airfoil  subject  to  turbulent  flow.  The  blue, 
solid  lines  represent  the  case  with  an  elastic  wall  while  the  black,  dashed  lines  are 
the  rigid  wall  case. 
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= 

Cv  PqVqvL 


(53) 


where  p0  is  the  free  stream  density,  t/0  is  the  free  stream  velocity,  L  is  the  length  of  the 

sensor,  v  is  the  square  root  of  the  autospectra  I  density  magnitude  of  the  "upwash"  compo¬ 
nent  of  the  turbulent  velocity,  and  GL>  is  the  square  root  of  the  autospectral  density  of  the 
unsteady  lift.  The  frequency  scale  was  non-dimensionalized  by  U0  and  the  integral  length 
scale  of  the  turbulence,  A.  defined  as  [8  j 

QO 

A  =  (54) 

o 

where  Rxx( t)  is  the  autocorrelation  function.  At  the  measurement  location  near  the  tip, 

A  =  7, 1  mm.  The  integral  length  scale  at  the  center  of  the  test  section  was  found  to  be  1 .2 
cm.  The  lift  coefficient  for  the  unsteady  lift  response  is  shown  in  Fig.  38  for  both  the  high¬ 
est  and  lowest  flow  speeds,  both  with  a  rigid  and  an  elastic  wall.  The  data  show  that  the 
elastic  wall  did  not  have  an  effect  on  the  unsteady  lift  response  of  the  airfoil.  The  only 
exception  was  at  a  non-dimensional  frequency  of  approximately  0.2,  where  the  unsteady 
lift  with  a  rigid  wall  showed  a  negative  peak  in  the  spectral  density.  However,  the  local 
minimum  was  at  a  constant  dimensional  frequency  and  therefore  did  not  change  with 
speed  which  suggests  that  this  result  was  not  a  physical  response  but  merely  an  artifact  of 
the  signal  processing,  likely  due  to  a  calibration  error.  Aside  from  that  one  notable  fre¬ 
quency.  the  data  in  Fig.  38  agree  remarkably  well.  Thus  the  unsteady  lift  can  be  considered 
to  be  independent  of  both  Reynolds  number  and  the  wall  condition  used. 
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6.4  Wall  Pressure  Spectrum 

The  wall  pressure  spectrum  was  measured  using  16  surface  mounted  microphones 


Figure  39:  Autospectral  density  of  the  two  microphones  located  nearest  to  the  cen 
terline  of  the  tunnel.  The  blue,  solid  curves  represent  the  case  with  the  airfoil 
present  while  the  black,  dashed  lines  are  the  case  with  no  airfoil  present. 


in  a  vertical  array,  embedded  in  a  rigid  wall.  The  microphones  were  placed  1 .9  cm  apart 


B 


Figure  40:  A.)  Magnitude  and  B.)  phase  of  the  cross  spectral  density  betw  een  the  two 
microphones  located  closest  to  the  centerline  of  the  test  section.  The  solid,  blue  lines 
represent  the  case  in  w  hich  the  airfoil  is  present  while  black,  dashed  lines  are  without 
the  airfoil. 
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Figure  41;  A.)  Magnitude  and  B.)  phase  of  the  cross  spectral  density  between  the  two 
microphones  located  furthest  from  the  centerline  of  the  test  section.  The  solid,  blue 
lines  represent  the  case  in  which  the  airfoil  is  present  while  black,  dashed  lines  are 
without  the  airfoil. 

with  8  microphones  on  each  side  of  the  centerline  of  the  test  section  as  shown  in  Fig.  37. 
Figure  40  shows  the  autospectral  density  of  the  wall  pressure,  normalized  by  the  dynamic 


pressure  (p{)IJ2q/2)  .  Four  spectra  are  shown  in  Fig.  40  which  represent  the  two  micro¬ 
phones  closest  to  the  centerline.  Measurements  with  and  without  the  airfoil  present  are 
shown.  The  results  indicate  that  the  presence  of  the  airfoil  had  a  small,  if  not  negligible 
effect  on  the  wall  pressure  spectra  for  frequencies  /A/  U{)  >  0.2 .  At  lower  frequencies  the 
airfoil  caused  a  substantial  increase  in  the  unsteady  pressure.  It  is  interesting  to  note  that 


both  configurations  indicate  an/  relationship  at  low  frequencies. 

Further  insight  into  the  effect  of  the  airfoil  on  the  wall  pressure  can  be  obtained 
from  the  cross  spectral  density  of  two  microphone  signals.  Figure  39  shows  the  magnitude 
and  phase  of  the  cross  spectral  density  betw  een  the  same  microphones  used  to  compute 
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the  auto  spectra  in  Fig.  40,  that  is,  0.95  cm  on  either  side  of  the  airfoil  leading  edge.  The 
magnitude  of  the  cross  spectrum  was  increased  by  as  much  as  a  factor  of  10  at  lower  fre¬ 
quencies,  with  negligible  increase  at  higher  frequencies.  The  phase  relationship,  shown  in 
Fig.  39b  indicates  that  for  low  frequencies  the  two  signals  were  in  phase  when  the  airfoil 
was  not  present,  and  n  out  of  phase  when  the  airfoil  was  present.  The  phase  matching 
without  the  airfoil  is  intuitive  given  that  the  microphones  were  placed  2.67 A  apart.  That 
is,  it  would  be  expected  that  the  turbulent  motions  that  contribute  to  the  spectra  at 
fA/UQ  «  i  have  a  length  scale  that  is  significantly  greater  than  the  microphone  spacing, 
and  thus  the  two  locations  are  subjected  to  essentially  the  same  low  frequency  unsteadi¬ 
ness.  The  switch  in  phase  to  -n  when  the  airfoil  is  present  further  suggests  that  the  pres¬ 
sure  field  is  dominated  by  the  turbulence-leading  edge  interaction.  Specifically,  both  the 
near  field  hydrodynamic  pressure  and  the  radiated  sound  pressure  will  be  dipole-like  due 
to  the  net  reaction  force  exerted  by  the  airfoil  onto  the  fluid. 

Figure  41  shows  the  cross  spectral  density  magnitude  and  phase  for  the  two  micro¬ 
phones  located  at  ±  14.3  cm.  These  were  the  furthest  away  from  the  centerline  of  the  test 
section.  These  show  the  same  general  behaviors  as  those  seen  in  Fig.  39,  however,  there 
are  two  distinct  differences.  First,  the  magnitude  of  the  correlation  between  these  two 
locations  is  a  factor  of  5  lower  than  that  for  the  two  microphones  closest  to  the  centerline. 
This  is  to  be  expected  because  the  spacing  of  these  two  microphones  is  much  greater  than 
the  integral  turbulent  length  scale.  Second,  since  the  hydrodynamic  field  at  these  locations 
was  not  likely  altered,  the  increase  in  the  wall  cross  spectral  density  w ith  the  airfoil 
present  can  be  attributed  to  acoustic  radiation  from  the  tip  of  the  airfoil.  The  phase  differ¬ 
ence  of  n  indicates  that  the  acoustic  radiation  from  the  airfoil  is  dipole  like  sound. 


University  of  Notre  Dame 


Center  for  Flow  Physics  and  Control 


Sound  and  Vibration  Measurements  Using  an  Air-analog  Model 


67 


6.5  Elastic  Wall  Vibration  Response 

The  elastic  wall  used  was  made  from  0.127  mm  thick  Aluminum  that  was  ten¬ 
sioned  isotropically  at  2 19N/m.  The  motion  of  the  elastic  wall  can  be  approximated  as  a 
solution  to  the  equation  for  waves  in  a  membrane  [  I  ], 

(V,  +  *2a2)ri  =  -P(r,  4>,  u)  (55) 

where  r\  is  the  displacement  of  the  membrane,  P  is  the  fluctuating  net  pressure  difference. 

k  =  co/c,  c  is  the  speed  of  sound,  and  a  =  c/  JTTa  ,  where  a  is  the  density  per  unit 
area,  and  T  is  the  tension  per  until  length.  The  effects  of  flexural  wave  motion,  aerody¬ 
namic  and  acoustic  damping,  and  fluid  loading  have  been  ignored.  For  the  case  of  a  round 
membrane  with  a  clamped  boarder,  the  solution  can  be  given  in  terms  of  the  eigenfunc¬ 
tions  for  the  free  vibration  of  the  membrane: 

cos  (mfy)Jm(2nvnr/ c  )  =  0  (56) 

sin 

where  Jm  is  the  Bessel  function  of  order  m,  r  is  the  radius  of  the  membrane,  cp  is  the  wave 
speed  of  the  membrane  and  vn  is  the  frequency  of  the  wave.  The  zeros  of  the  Bessel  func¬ 
tion  can  be  found  to  provide  the  natural  frequencies  of  oscillation  allowed  in  the  mem* 
brane.  For  the  experimental  conditions  given  here,  the  natural  frequencies  are  shown  in 
Table  I ,  The  structural  response  of  the  membrane  at  two  locations,  both  with  and  without 
the  airfoil,  is  shown  as  an  auto  spectral  density  of  the  membrane  velocity  in  Figs.  42  A  and 
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42 B.  The  physical  locations  chosen  for  these  figures  were  somewhat  arbitrary,  but  the  data 
are  indicative  of  observations  at  other  locations  as  well. 


Table  1:  Natural  frequencies  of  vibration  for  a  clamped,  circular  membrane. 


n/m 

0 

1 

2 

3 

4 

5 

0 

0.014 

0.022 

0.029 

0.036 

0.043 

0.049 

1 

0.031 

0.039 

0.047 

0.055 

0.062 

0.069 

2 

0.049 

0.057 

0.065 

0.073 

0.081 

0.088 

3 

0.066 

0.075 

0.083 

0.091 

0.099 

0.107 

4 

0.084 

0.093 

0.101 

0.109 

0.117 

0.125 

The  addition  of  the  airfoil  was  observed  to  create  a  substantial  increase  in  the 
vibration  amplitudes  shown  in  Fig.  42  A.  The  frequency  range  and  amplitude  of  the 
increase  were  generally  consistent  with  the  increased  auto  and  cross  spectral  measure¬ 
ments  of  the  wall  pressure  shown  in  Figs.  39  and  4 1 .  Specifically,  the  rough  increase  in  the 
wall  pressure  by  a  factor  of  10  in  the  range  0.01  </A/(/0  <  0. 1 5  can  be  observed  in  Fig. 
42A  as  well.  However,  it  is  clear  that  the  increased  phase  organization  of  the  wall  pressure 


Figure  42:  Autospectral  density  of  the  vibrational  response  of  the  membrane  at  A.) 
a  radius  of  14.7  cm  and  an  angle  of -105  deg  and  B.)  a  radius  of  12.19  cm  and  an 
angle  of  90  deg.  The  solid,  blue  curve  is  with  the  airfoil  present  while  the  black, 
dashed  curve  is  in  the  absence  of  the  airfoil.  +  indicates  position  on  membrane. 
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due  to  the  airfoil  was  causing  some  of  the  modes  to  be  excited  at  significantly  higher  lev¬ 
els.  For  example,  the  fundamental  membrane  mode  (m=n- 0)  was  increased  by  a  factor  of 
nearly  100.  At  higher  frequencies  (/A/f/0>  0.1 )  the  vibration  response  was  nearly  iden¬ 
tical  for  the  two  test  cases. 

The  data  shown  in  Fig.  8  were  acquired  under  the  same  conditions  as  those  in  Fig. 
42A.  but  show  a  number  of  notable  differences.  First,  the  broad  increase  in  amplitudes  at 
lower  frequencies  observed  in  Figure  42A  is  not  present  in  Fig.  42B.  In  addition,  modes 
that  were  not  excited  with  the  homogeneous  turbulent  forcing,  for  example  [0, 1 )  and  10,2], 
were  found  to  be  excited  by  the  presence  of  the  airfoil.  At  higher  frequency  the  vibration 
amplitude  was  found  to  by  increased  by  a  factor  of  5-10  in  the  range 
0.07  </A/f/0  <  0.25 .  Given  that  the  magnitude  of  the  wall  pressure  fluctuations  were 

found  to  be  unaffected  by  the  airfoil  at  these  frequencies,  it  is  reasonable  to  expect  that  the 
relative  phase  of  the  wall  pressure  led  to  this  increase. 


Figure  43:  Autospcctral  density  of  a  single  far  field  microphone  plotted  on  a  Deci¬ 
bel  scale  with  a  bin  width  of  6.9x10-4  vs  non-dimensional  frequency.  The  solid,  blue 
line  is  the  case  with  the  airfoil  while  the  black,  dashed  line  is  without  the  airfoil 
present. 
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6.6  Acoustic  Noise  Response 

The  acoustic  response  of  the  elastic  wall  to  complex  turbulent  forcing  can  be  rep¬ 
resented  as  the  convolution  of  the  membrane  motion  with  the  structural  Green's  func¬ 
tion.!  1]  The  near  field  is  dominated  by  evanescent  waves  which  cancel  in  the  far  field.  The 
modes  that  do  not  cancel  are  scattered  off  of  the  boundaries  of  the  membrane  or  other  sur¬ 
face  discontinuities.  These  non-cancelling  modes  will  radiate  to  the  far  field  as  acoustic 
waves. 

Data  were  taken  with  a  single  far  field  (2  m)  microphone  both  with  and  without  the 
airfoil  present,  figure  43  shows  the  far  field  acoustic  spectra  for  both  cases.  Over  a  range 
of  0.04  <  fA/Uy  <0.1,  the  airfoil  increased  the  far  field  acoustic  noise  by  approximately 

2  dB.  An  increase  of  this  nature  is  to  be  expected  given  the  response  of  both  the  wall  pres¬ 
sure  spectra  and  elastic  wall  vibration  response  to  the  presence  of  the  airfoil.  I  lowever,  the 
range  of  frequencies  over  which  there  is  a  substantial  increase  is  smaller  than  (hat 
observed  for  the  wall  pressure  and  elastic  wall  vibration  spectra.  There  was  also  a  large 
increase  of  approximately  3  dl3  at  a  frequency  of  0.22,  which  is  very  close  the  natural 
shedding  frequency  of  the  turbulence  grids. 

6.7  Conclusions 

The  sound  and  vibration  of  an  elastic  wall  produced  by  an  airfoil  tip.  in  a  turbulent 
How,  near  the  wall  was  investigated.  An  experiment  was  conducted  using  two  sets  of 
boundary  conditions,  the  first  being  a  canonical  setup  in  which  the  membrane  is  excited  by 
turbulent  flow  alone  and  the  second  adding  a  flat  plate  airfoil  to  the  flow  with  a  narrow  tip 
gap  to  the  elastic  wall.  The  unsteady  lift  response  of  the  airfoil  was  measured  using 
unsteady  force  transducers.  It  was  determined  that  the  elastic  wall  had  a  negligible  effect 
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on  the  unsteady  lift  response  when  compared  to  that  with  a  rigid  wall  in  place.  The  pres¬ 
ence  of  the  airfoil  had  a  significant  effect  on  the  wall  pressure  spectrum  over  a  range  of 
0.0 1  <  f A/ Uq  <  0.1 5 .  The  effect  was  also  observed  in  the  vibrational  response  of  the  elas¬ 
tic  wall.  The  increased  structural  vibration  of  the  membrane  corresponded  to  a  1-2  dB 
increase  in  the  far  field  radiated  acoustics.  This  was  a  significant  increase  in  the  far  field 
sound  generated  by  an  elastic  wall  alone. 
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7.0  RECIPROCITY  RESULTS  FOR  RIGID  AND  ELASTIC  DUCTS 

The  study  of  sound  radiated  from  aeroacoustic  and  structural  acoustic  systems  requires  an  under¬ 
standing  of  both  the  acoustic  sources  as  well  as  the  transfer  of  acoustic  energy  from  the  source 
location  to  a  point  of  interest  (e.g.,  the  far  field).  The  experimental  use  of  reciprocity  principles  to 
measure  the  acoustic  Green's  function  was  described  in  Section  4.  This  section  presents  results 
from  the  reciprocal  measurement  of  the  acoustic  Green’s  function  and  its  gradient  for  two  geome¬ 
tries  of  interest.  The  first  was  a  rigid  PVC  cylindrical  duct  with  8. 1  inch  inner  diameter  and  24 
inch  length.  The  second  was  constructed  from  three  sections  of  the  thin  nickel  shell  described  in 
Section  3.  The  end  conditions  for  both  the  rigid  and  elastic  ducts  were  partially  flanged  as  shown 
in  Figure  44.  The  internal  coordinate  system  shown  describes  the  (r,9,z)  location  of  the  virtual 
acoustic  source  in  the  duct.  The  far  field,  or  “receiver”  location  is  shown  in  Figure  45. 

The  purpose  of  the  measurements  was  to  investigate  the  dependence  of  the  far  field  sound  on  the 
source  and  far  field  locations,  and  to  compare  these  results  for  the  rigid  and  clastic  ducts.  A 
related  objective  was  to  determine  if  the  sensitivity  of  the  far  field  sound  to  the  source  location  is 
different  for  the  elastic  duct  as  compared  to  the  rigid  boundary  conditions  with  nominally  the 
same  geometry.  This  issue  is  important  for  the  definition  of  blocked,  or  net  equivalent  acoustic 
sources.  Specifically,  if  a  distribution  of  sources  (such  as  a  ducted  rotor)  is  modeled  as  a  single 
point  source  in  a  rigid-walled  system,  significant  errors  may  be  encountered  when  the  actual 
sources  are  placed  in  an  system  where  the  elastic  walls  may  change  the  spatial  sensitivity  of  the 
local  Green's  function. 

Measurements  were  conducted  as  described  in  Section  4  to  obtain  the  Green's  function,  G,  as  well 
as  its  gradient,  GD.  These  functions  represent  the  acoustic  transfer  functions  for  monopole  and 
dipole  sources,  respectively.  For  ease  of  presentation,  the  results  presented  in  the  following  sub¬ 
sections  will  be  given  as  the  difference  between  the  measured  value  and  the  magnitude  that  would 
be  obtained  in  a  free  space  environment.  Specifically. 

G'(V,)S  lOloglGt*^)]-  lOlogtG^H)]  (57) 

where  GfS(|H)  is  the  free  space  Green's  function  evaluated  at  the  radial  location  |r|  =  |x0-.vv| . 

Similarly,  the  gradient  function  difference  is  defined  as 

Gb(x0,xsy=  lO|og([GD(x0,^)])-IO]og[GD.FS(|r|)]  (58) 

w  here  GD  is  defined  in  equation  7,  and  the  free  space  value,  GI}.[tS  is  defined  in  equation  23. 
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Figure  45  Far-field  "receiver"  coordinate  system 
7.1  Far  field  angular  dependence  of  Green’s  function 

This  section  presents  results  that  describe  the  angular  dependence  of  the  G  and  G^.  The  source 
location,  xs,  was  placed  at  the  origin  of  the  coordinate  system  shown  in  Figure  44.  The  tar  field 
location  was  measured  2m  from  the  duct  center,  with  the  angular  coordinate.  <f>.  treated  as  the 
independent  variable  of  interest. 

Contours  of  O'  and  Gp  are  shown  in  Figures  46  through  49  for  both  the  rigid  and  elastic  ducts.  In 

these  "polar"  plots  the  radial  direction  represents  the  frequency,  and  the  angular  position  repre¬ 
sents  the  far  field  angle.  <|).  For  example,  the  values  shown  near  the  horizontal  axis  of  the  polar 
plot  represents  the  values  "on  axis”  with  the  duct.  The  vertical  direction  represents  a  far  field  loca¬ 
tion  that  is  perpendicular  to  the  duct  axis.  The  frequency  values  are  shown  in  both  axes  scaled  by 
the  fundamental  organ  pipe  frequency  (eqn  24)  on  the  horizontal  axis,  and  normalized  by  the  duct 
radius  (a)  on  the  vertical  axis.  Also,  the  frequencies  of  higher  order  mode  cut-on  are  noted  by  the 
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ordered  pairs  in  the  figure  where  the  first  integer  represents  the  radial  mode  order,  and  the  second 
represents  the  circumferential  mode  order,  for  the  G'D  measurements,  the  dipole  axis  was  aligned 
with  the  duct  axis. 

The  contours  of  G'  shown  in  Figure  46  illustrates  a  number  of  distinctive  features  in  the  radiated 
sound  field.  First,  the  first  organ  pipe  mode  is  observed  to  resonate  which  leads  to  a  lOdB  increase 
in  far-field  sound  at  all  angles  with  respect  to  the  duct.  The  second  pipe  mode  is  not  excited  since 
the  axial  center  is  an  anti-node  for  this  frequency.  The  third  pipe  mode  is  again  excited,  although 
two  distinctive  peaks  can  be  observed  in  the  azimuthal  direction.  The  first  is  on  axis,  and  the  sec¬ 
ond  is  approximately  60  degrees  from  axis.  The  radiation  at  higher  frequency  is  observed  to  be 
primarily  on  axis,  with  some  increased  levels  between  20  and  60  degrees.  The  cut-on  resonances 
are  notable,  and  show  several  local  maxima  in  the  azimuthal  direction. 

The  angular  dependence  of  %  for  the  rigid  duct  is  shown  in  Figure  47.  These  data  show  similar 
trends  compared  to  the  Green's  function  described  above.  A  notable  difference  is  the  excitation  of 
the  even,  rather  than  the  odd  ordered  pipe  modes.  This  is  expected  since  the  anti-nodes  of  the  gra¬ 
dient  function  will  be  located  at  the  nodes  ofthe  original  function.  It  is  also  noted  that  the  cut-on 
frequencies  are  not  amplified  (compared  to  the  free  space  values). 

Contour  plots  of  G'  and  Gq  for  the  elastic  duct  are  shown  in  Figures  48  and  49  respectively,  it  is 
first  noted  that  the  general  character  ofthe  sound  radiation  with  the  elastic  duct  is  similar  to  the 
rigid  case.  Specifically,  local  maxima  in  the  rigid  duct  functions  generally  correspond  to  local 
maxima  for  the  elastic  duct.  I  lowever,  it  is  observed  that  generally  increased  radiation  levels  are 
present  with  the  elastic  duct  sections.  This  is  particularly  evident  at  higher  frequencies  where  the 
higher-order  duct  acoustic  modes  are  interacting  with  the  vibrating  duct. 
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Figure  46  Contours  of  G*  showing  the  angular  and  frequency  dependence  for  the  rigid  ducL 
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Figure  48  Contours  of  G'  showing  the  angular  and  frequency  dependence  for  the  elastic  duct. 


r 

Figure  49  Contours  of  %  showing  the  angular  and  frequency  dependence  for  the  rigid  duct. 
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7.2  Effect  of  Axial  Location  of  Source  on  Radiated  Sound 

The  effects  of  the  axial  location  of  the  sound  source  are  presented  in  Figures  50  through  53.  The 
contour  plots  are  shown  with  the  frequency  on  the  abscissa,  again  normalized  by  both  the  first 
organ  pipe  frequency  as  well  as  the  duct  radius  on  the  lower  and  upper  axes,  respectively.  The 
ordinate  of  the  plots  shows  the  axial  location  normalized  by  the  duct  length  (-0.5<x/L<0.5). 

The  magnitude  of  G’  for  the  rigid  case  shown  in  Figure  50  indicates  a  number  of  distinctive  and 
interesting  features  to  the  radiated  sound  field.  First,  the  nodal  pattern  of  the  first  organ  pipe  mode 
can  be  observed  by  the  strong  amplification  at  the  duct  center  compared  to  the  duct  ends  at  /*=l . 
Similarly,  at_/*=2  the  amplifications  are  observed  at  axial  locations  of  -0.25  and  0.25  correspond¬ 
ing  to  the  nodes  of  the  standing  wave.  Note  that  the  dependence  on  z/L  is  not  symmetric  about  the 
duct  center  since  the  far  field  location  is  placed  on  axis  on  the  z>0  side  of  the  duct.  The  magnitude 
of  Gf}  shown  in  Figure  5 1  for  the  rigid  duct  shows  similar  results,  although  the  organ  pipe  modes 
are  excited  at  the  anti-nodes  of  the  standing  wave  as  expected. 

The  radiated  sound  exhibits  dramatic  changes  in  character  near  both  the  (0. 1 )  and  (0.2)  cut  on  fre¬ 
quencies.  It  is  clear  that  the  magnitudes  of  G'  and  Gfo  at  high  frequency  is  both  considerably 
larger  than  the  free  space  Green’s  function,  and  is  also  highly  dependent  on  the  location  of  the 
source.  For  example,  at  /*=I3,  G'  varies  by  nearly  20dB  from  z/L=0.5  to  0.3. 

The  results  for  the  elastic  duct  generally  exhibit  increased  magnitudes  at  all  frequencies  and  axial 
locations.  It  is  also  noted  that  there  are  fewer  locations  at  which  the  Green's  function  and  its  gra¬ 
dient  have  a  lower  magnitude  than  the  free  space  values.  Note  that  the  structural  free  vibration  fre¬ 
quencies  of  the  shell  are  all  about  200Hz,  which  corresponds  to/*=0.7.  Issues  related  to  the 
structural  vibration  will  be  described  in  the  next  section. 
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Figure  50  Axial  dependence  of  G'  for  rigid  duct 
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Figure  5 1  Axial  dependence  of  G rD  for  rigid  duct 
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Figure  52  Axial  dependence  of  G '  for  elastic  duct 
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Figure  53  Axial  dependence  of  G 'D  for  elastic  duct 
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7.3  Effect  of  Radial  Location  of  Source  -  On  and  OfT  Axis 

The  radial  position  of  an  axial  source  is  an  important  parameter  in  predicting  the  radiated  sound 
from  a  ducted  rotor  system.  Specifically,  a  ducted  rotor  will  consist  of  a  radial  distribution  of 
sources  that  may  lead  to  varied  levels  of  plane  wave  and  higher  order  mode  excitation  within  the 
duct.  If  the  distributed  acoustic  sources  are  modeled  as  a  single,  net  source,  then  accurate  predic¬ 
tions  of  radiated  sound  can  only  be  made  if  an  appropriate  weighted  average  of  the  acoustic 
Green’s  function  is  made. 

The  results  showing  G'  and  Gp  for  the  rigid  duct  with  tar  field  location  on  axis  with  the  duct  are 
shown  in  Figures  54  and  55.  respectively.  The  axial  location  of  the  source  was  z/L=0.  The  same  is 
plotted  for  the  elastic  duct  in  Figures  56  and  57.  These  four  plots  are  repeated  with  the  far  Held 
location  at  and  angle  of  45  degrees  from  the  duct  axis  in  Figures  58  through  61.  In  all  figures  the 
horizontal  axes  represent  the  frequency  normalized  by  either  the  first  organ  pipe  frequency  or  as 
ka  as  before.  The  ordinate  represents  the  radial  position  normalized  by  the  duct  radius. 

The  rigid  duct  exhibits  nearly  uniform  magnitude  of  both  G'  and  G p  in  the  radial  direction  at 
lower  frequencies  where  the  plane  wave  mode  is  dominant.  As  expected  from  the  axial  location  of 
the  source,  of  G'  is  dominated  by  the  odd  ordered  organ  pipe  modes  at  all  radial  locations, 
whereas  Gp  is  dominated  by  the  even  ordered  modes.  At  higher  frequencies  the  higher  order 
modes  lead  to  more  substantial  radial  dependence  of  both  G'  and  Gp.  For  example,  the  (0. 1 ), 

( I .  I ).  and  (0,2)  modes  appear  to  be  more  strongly  excited  at  the  duct  center  (near  r/a= 0).  The 
modes  (3,0)  and  (3, 1 )  are  considerably  stronger  near  the  duct  wall  (r/a=  1). 

A  number  of  interesting  features  are  observed  when  comparing  the  results  for  the  elastic  shell. 
First,  the  lower  frequencies  are  observed  to  exhibit  radial  dependence  that  was  not  observed  for 
the  rigid  duct.  The  on  axis  data  (Figures  56  and  57)  indicate  an  increased  sensitivity  to  the  radial 
position  at  higher  frequency,  with  particularly  large  magnitudes  observed  near  the  duct  centerline. 
At  45  degrees  off  axis,  (Figures  60  and  61 ),  the  elastic  shell  exhibits  the  largest  values  at  the  outer 
radial  locations  at  high  frequency. 


University  of  Notre  Dame 


Center  for  Flow  Physics  and  Control 


Sound  and  Vibration  Measurements  Using  an  Air- analog  Model 


81 


Figure  54  Radial  dependence  of  G'  (on  axis)  for  rigid  duct 
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Figure  55  Radial  dependence  of  G*D  (on  axis)  for  rigid  duct 
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Figure  56  Radial  dependence  of  G'  (on  axis)  for  elastic  duct 


Figure  57  Radial  dependence  of  G’D  (on  axis)  for  elastic  duct 
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Figure  58  Radial  dependence  of  G'  (45  degrees  from  axis)  for  rigid  duct 
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Figure  59  Radial  dependence  of  G'D  (45  degrees  from  axis)  for  rigid  duct 
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Figure  60  Radial  dependence  of  G'  (45  degrees  from  axis)  for  elastic  duct 


Figure  61  Radial  dependence  of  G„  (45  degrees  from  axis)  for  elastic  duct 
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8.0  PROPELLER  -  ELASTIC  DUCT  INTERACTIONS 

8.1  Introduction 

The  original  objectives  of  this  research  involved  the  development  of  the  experimental  tools  neces¬ 
sary  to  investigate  the  interaction  between  the  acoustic  sources  of  a  ducted  rotor,  the  convecting 
hydrodynamic  pressure  field,  and  the  shell  vibrations.  Consultation  with  engineers  from  NSWC. 
ONR,  and  Electric  Boat  provide  significant  feedback  and  technical  guidance  for  this  part  of  the 
research.  One  of  the  primary  objectives  that  was  added  during  the  last  phase  of  the  project 
involved  the  investigation  of  inlet  distortion  effects  on  the  vibro-acoustic  response  of  a  ducted 
rotor.  This  results  presented  in  this  section  are  intended  to  be  a  status  update  for  this  work,  which 
is  continuing  under  ONR  Grant  #NG00 1 4-07-1  -0570. 

The  general  boundary  conditions  and  experimental  setup  for  the  propeller-duct  interaction  experi¬ 
ments  were  described  in  Section  2.  The  duct  length  was  3 1  inches,  with  a  corresponding  first 
organ  pipe  mode/*=i  89Hz.  The  significant  change  of  relevance  to  the  present  work  is  the  use  of 
varied  inlet  section  design  such  that  variable  inlet  total  pressure  as  well  as  the  overall  level  and 
spectral  content  of  the  turbulent  fluid  motions  were  obtained.  These  inlet  conditions  were  pro¬ 
vided  in  both  an  axisymmetric  as  well  as  a  non-symmetric  configuration  as  described  in  the  fol¬ 
lowing  subsections.  These  experiments  provide  a  method  to  directly  access  the  impact  of  inlet 
distortion  on  the  vibration  response  of  the  shell. 

8.2  Rotor  Inflow  Boundary  Conditions 

The  properties  of  the  mean  flow  as  well  as  the  turbulence  were  documented  in  detail  for  several 
inlet  duct  section  geometries.  The  measurements  were  acquired  using  standard  hot-wire  anemom- 
etry  techniques.  Specifically,  a  single  hot-wire  probe  was  calibrated  and  traversed  radially  in  the 
duct  at  the  rotor  inlet  plane.  Airflow  was  provided  by  an  auxiliary  blower  for  these  experiments. 
The  mean  velocity  and  RMS  of  the  fluctuating  velocity  acquired  for  the  baseline,  smooth  inlet 
geometry  is  shown  in  figure  62.  Data  are  shown  for  several  streamwise  locations  with  respect  to 
the  duct  inlet.  The  data  labeled  L0  were  acquired  at  the  inlet  of  the  duct,  and  exhibit  a  typical 
inviscid  ‘‘overshoot"  due  to  the  strong  streamline  curvature  that  exists  at  the  inlet.  The  location 
“L2"  (2  inches  from  the  inlet)  shows  a  more  standard  velocity  profile.  The  boundary  layer  thick¬ 
ness  is  observed  to  increase  monotonically  with  increased  distance  downstream  as  expected.  The 
corresponding  RMS  values  are  consistent  with  typical  boundary  layer  measurements. 
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The  autospectral  densities  from  several  locations  in  the  boundary  layer  are  shown  in  Figure  63, 
These  data  are  important  because  they  form  the  input  excitation  for  the  aeroacoustic  source  terms. 
It  should  be  observed  that  the  spectral  levels  continue  to  increase  at  low  frequency.  That  is,  there 
is  no  “leveling  off’  of  the  spectra  at  low  frequency  that  would  normally  be  expected  in  free  shear 
turbulence.  The  reason  for  this  is  well  understood  in  the  boundary  layer  turbulence  community, 
and  is  related  to  the  long  stream  wise  structures  that  exist  in  both  the  very  near  wall  region  as  well 
as  the  outer  regions  of  the  boundary  layer.  This  is  an  important  aspect  for  modeling  the  hydroa¬ 
coustic  and  structural  acoustic  problem.  Specifically,  the  majority  of  the  unsteady  fluid  motions 
that  lead  to  structural  vibrations  can  be  considered  to  be  strongly  anisotropic  in  nature.  The  higher 
frequency  (small  scales)  that  are  more  closely  model  by  isotropic  turbulence  theory  are  consider¬ 
ably  smaller  in  magnitude  and  thus  to  not  contribute  significantly  to  the  acoustic  signature  of  the 
ducted  rotor  system. 


Figure  62  Mean  velocity  and  RMS  of  fluctuating  velocity  acquired  at  various  streamwise  loca¬ 
tions  (in  inches  from  the  inlet)  for  the  smooth  wall  inlet  geometry. 
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Figure  63  Sample  velocity  spectral  results  from  smooth  inlet  boundary  layer  at  a 
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A  second  rotor  inlet  geometry  was  designed,  machined,  and  tested  in  the  ducted  rotor  facility.  A 
solid  model  image  of  the  geometry  is  shown  in  Figure  64.  The  interior  wall  of  this  section  is  con¬ 
figured  as  the  standard  "NASA  hump”.  The  profile  of  the  NASA  hump  was  designed  specifically 
to  have  a  fixed  separation  point  on  a  convex  surface.  The  advantage  of  this  geometry  for  the 
present  application  is  that  it  provides  a  well  defined,  large  scale  turbulent  approach  flow  to  the 
rotor  w  ithout  generating  additional  noise.  Thai  is,  since  the  hump  profile  does  not  contain  sharp 
comers  it  will  not  radiate  sound  efficiently,  and  thus  will  not  add  additional  sources  to  be  consid¬ 
ered  in  the  hydro/structural  problem. 

The  velocity  results  from  the  symmetric  NASA  hump  are  shown  in  Figures  65  and  66.  The  mean 
and  RMS  of  velocity  describe  a  typical  "wake-like”  behavior  that  is  followed  by  a  redeveloping 
boundary  layer  at  the  downstream  locations.  The  spectral  results  indicate  a  prominent  shedding 
frequency  that  decays  in  the  stream  wise  direction.  The  furthest  location  measured  (L8)  indicates 
that  the  turbulence  has  "relaxed"  to  a  state  that  can  be  roughly  approximated  as  isotropic  turbu¬ 
lence.  That  is.  three  distinct  regions  can  be  identified  in  the  spectral  plot:  a  flat  region  at  low  fre¬ 
quency,  a  -5/3  slope  region  characteristic  of  the  inertial  subrange,  followed  by  an  exponential  roll¬ 
off  at  the  highest  frequencies. 
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Figure  65  Mean  and  RMS  of  velocity  in  the  wake  of  the  NASA  hump 
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Figure  66  Sample  velocity  spectra  obtained  in  the  wake  of  the  NASA  hump. 

The  third  and  final  geometry  constructed  for  this  study  is  shown  in  Figure  67.  The  models  is  con¬ 
figured  to  generate  a  circumferential  distortion  to  the  rotor.  Specifically,  the  height  of  the  NASA 
hump  described  in  the  previous  geometry  was  varied  around  the  circumference  from  zero  at  one 
location  (say,  0  degrees)  to  a  maximum  of  0.375  inches  (same  as  above)  at  180  degrees.  This 
resulted  in  a  circumferential  varied  inlet  velocity  and  inlet  turbulent  level  as  shown  in  Figure  68. 
The  corresponding  spectral  results  are  shown  in  Figure  69. 


Figure  67  Model  of  the  non-symmetric  inlet  geometry 
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Figure  68  Mean  and  RMS  contours  downstream  of  the  non-symmetric  inlet. 
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Figure  69  Sample  velocity  spectra  downstream  of  the  non-symmetric  hump 


University  of  Notre  Dame 


Center  for  Flow  Physics  and  Control 


Sound  and  Vibration  Measurements  Using  an  Air -analog  Model 


90 


8.3  In-vacuo  Analysis  of  Elastic  Shell 

The  natural  frequencies  of  the  elastic  shell  can  be  determined  from  analysis  of  the  Donnell  equa¬ 
tions.  By  assuming  the  number  of  circumferential  harmonics  is  large  compared  to  the  axial  wave- 
number  (Yu's  assumption),  the  natural  frequencies  are  given  by  equation  2.88  of  Leissa  (1973): 

fi2  =  *«*  +  (!  -u)2(^)\  (59) 

The  nondimensional  frequency,  fi  is  given  by: 

n2  =  e.U_.y.)fty.  (60) 

where  E  is  the  elastic  modulus  of  the  shell,  and  to  is  the  circular  frequency  (2n/).  The  axial  wave¬ 
number  is  given  by  the  parameter  Xm.  defined  as: 


where  em  represent  the  roots  of  the  clamped-c  lamped  shell  equation: 


(61) 


given  as 


cos  (s,„)  cosh  (sm)-  I  -  0 

Bm  =  l.506n,  2.500rc,  3.500n.  etc... 


(62) 


The  results  for  the  parameters  given  for  the  thin  nickel  shell  are  shown  in  Figure  70.  The  modal 
density  is  shown  in  Figure  71. 

The  natural  frequencies  of  the  shell  were  determined  experimentally  by  exciting  the  shell  using  a 
speaker  in  close  proximity  to  the  shell.  The  cross  spectral  density  between  the  shell  vibration  sig¬ 
nal  and  the  excitation  signal  provided  the  amplitude  and  phase  information  needed  to  determine 
the  shell  mode  order  by  inspection.  Contour  plots  of  the  vibration  modes  are  shown  in  Figure  72. 
Each  image  shows  120  degrees  of  the  shell  that  represents  the  practical  limit  of  the  projected  view 
of  the  shell  from  a  stationary  laser  source. 

It  is  clear  from  inspection  of  Figure  70  that  the  true  frequencies  are  lower  than  the  values  pre¬ 
dicted  by  the  in-vacuo  shell  theory.  It  is  presumed  that  this  is  due  to  the  fluid  loading  effects  on 
the  shell.  Specifically,  at  300Flz,  the  fluid  loading  factor  (pe)/(ps/Ko)=0.35.  It  can  also  be 

noticed  that  two  of  the  circumferential  modes  ( 1 5  and  16)  have  "locked  on”  to  the  organ  pipe  fre¬ 
quency  of  189  Flz.  The  fluid  loading  factor  at  this  frequency  is  0.58,  which  suggests  that  the  fluid 
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loading  stress,  the  extensional  (membrane)  and  flexural  stresses  are  all  important  in  determining 
the  modal  frequencies. 

As  a  final  note,  the  axial  profile  of  the  mode  shapes  of  the  in-vacuo  shell  vibration  is  of  interest 
for  comparison  with  the  experimental  data.  The  radial  motion  is  given  by  (Leissa,  1973  eqn  2.86) 

w  =  2c[sinh(e)  -  sine - ; — - 1  .  (63) 

L  coshe-cosej 

[(sinhe  -  sine)(coshX.s-  cosXs)-(coshe  -  cose){sinhX..s-  sin^.v)] 

'['he  first  three  mode  shapes  are  shown  in  Figure  73. 


Figure  70  Theoretical  (in-vacuo)  shell  frequencies  (connected  for  reference)  and  mea¬ 
sured  values  (open  symbols). 
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Figure  71  Modal  density  for  the  elastic  shell 


Frame  003  !  15  Aug 2007  !  InsJ.  Vais®  24B.1H* MFramfeflOl  15Aug2007  IrwL  Vats @  184.4m 


University  of  Notre  Dame 


Center  for  Flow  Physics  and  Control 


Sound  and  Vibration  Measurements  Using  an  Air-analog  Model 


93 


Figure  73  Axial  profile  of  the  eigenmodes  of  a  clamped -clamped  shell 
8,4  Effects  of  Inlet  Distortion 

The  main  objectives  in  the  final  phase  of  this  research  involved  the  measurement  of  shell  vibra¬ 
tion  amplitudes  while  the  rotor  ingests  a  distorted  inlet  Mow.  Under  these  condition,  non-uniform 
convecting  hydrodynamic  pressures  are  present,  as  well  as  a  circumferential  varying  acoustic 
source  magnitude. 

The  experiments  were  conducted  as  follows.  First,  the  distorted  inlet  geometry  described  above 
was  installed  at  the  inlet  to  the  ducted  rotor  facility.  The  scanning  LDV  and  single  point  LDV 
were  then  used  to  obtain  the  two  point  correlation  measurements  of  the  duct.  The  non-symmetric 
portion  of  the  inlet  was  then  "clocked"  at  90  degree  intervals  with  respect  to  the  laser  location  and 
the  two  point  scan  was  repeated.  Preliminary  results  from  a  single  point  on  the  duct  are  shown  in 
Figures  74  and  75.  From  a  broad  view,  the  overall  amplitude  of  the  vibration  is  not  significantly 
different  at  the  four  locations  around  the  duct.  However,  some  clear  differences  are  noted.  First, 
the  low  frequencies  show  a  broadband  excitation  that  is  1  OdB  higher  for  two  of  the  neighboring 
locations.  Additionally,  the  closer  view  of  the  lowest  frequency  modes  in  Figure  75  suggest  that 
substantial  differences  do  exist  between  the  four  locations.  It  appears  that  several  of  the  modes 
have  either  shifted  frequency,  or  are  significantly  less  excited  at  certain  locations. 
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The  current  focus  is  being  placed  on  continued  analysis  of  the  acquired  data,  For  example,  the 
eigenmodes  of  the  full  two-point  correlation  matrix  will  allow  the  amplitudes  of  each  azimuthal 
mode  order  to  be  compared  for  the  various  conditions.  Additionally,  the  dependence  of  the  struc 
tural  vibration  on  rotor  speed  is  currently  being  examined. 


Figure  74  Autospectrai  densities  of  duct  vibration  at  4  azimuthal  locations. 
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Figure  75  Detailed  view  of  vibration  spectral  amplitudes  at  two  circumferential  loca¬ 
tions  of  the  duct. 


University  of  Notre  Dame 


Center  for  Flow  Physics  and  Control 


Sound  and  Vibration  Measurements  Using  an  A  ir -analog  Model 


95 


9.0  DISCUSSION  AND  CONCLUSIONS 

In  addition  to  the  detailed  conclusions  provided  in  each  of  the  sections  of  this  proposal,  it  is  useful 
to  provide  some  overall  comentary  and  conclusions  related  to  the  successes  of  the  research  pro¬ 
gram  as  well  as  issues  related  to  the  ongoing  work. 

The  initial  objective  of  this  research  was  to  develop  the  experimental  tools  required  for  conduct¬ 
ing  a  broad  range  of  studies  related  to  structural  acoustics.  More  specifically,  the  use  of  air  as  a 
working  fluid  for  investigating  fluid-loaded  cylindrical  shells  with  complex  hydrodynamic  and 
acoustic  forcing  was  considered.  The  challenge  in  this  regard  was  to  create  high  quality  test  spec¬ 
imens  that  were  very  thin,  thus  increasing  the  importance  of  fluid  loading.  Handeling  of  the  thin 
materials,  applying  uniform  and  measured  static  tension,  and  fixing  uniform  boundary  conditions 
were  all  important  considerations  for  the  developement  of  the  experimental  apparatus. 

The  facilities  and  instrumentation  already  avaialble  included  the  Anechoic  Wind  Tunnel,  an  80 
microphone  phased  array,  and  the  initial  design  of  the  rigid-walled  ducted  rotor  system.  During 
this  project,  two  laser  Doppler  vibrometry  systems  as  well  as  an  array  of  miniature  microphones 
were  added.  Data  acquisition  and  processing  systems  were  also  upgraded. 

As  described  in  this  report,  two  new  experimental  test  con figrat ions  were  designed  and  manufac¬ 
tured.  The  first  was  a  flat  plate  arrangment  in  which  various  materials  could  be  stretched  over  a 
template  with  an  arbitrary  opening.  This  allowed  the  study  of  how  the  fluid  loading  effects  the 
natural  frequencies  of  the  membrane. The  geometry  also  provided  a  uinquc  oportunity  to  investi¬ 
gate  the  effects  of  complex  hydrodynamic  forcing.  The  capabilities  developed  in  these  experi¬ 
ments  presents  a  unique  oportunity  for  studies  in  structural  acoustics  and  fluid-structure 
iterations.  Specifically,  the  effects  of  fluid  loading,  complex  and  non-homogeneous  convening 
pressures,  and  acoustic  forcing  can  all  be  considered  in  terms  of  modal  frequencies,  mode  shape 
functions,  and  the  net  radiated  sound  from  the  induced  structural  motion. 

The  second  new  facility  that  was  designed  and  constructed  under  this  project  was  the  nickel  shell 
used  in  the  ducted  rotor  rig.  Although  the  measurement  objectives  and  capabilities  are  similar  to 
the  flat  plate/membrane  appratus,  a  number  of  additional  and  complex  scientific  objectives  can  be 
addressed  due  to  the  complex  nature  of  the  shell  vibrations.  Also,  a  number  of  practical  problems 
can  be  addressed  such  as  the  the  interaction  with  swirling  rotor  w'akes,  the  fluid  loading  effects, 
etc.  The  present  work  began  a  new  project  where  these  features  are  examined  using  inflow  bound- 
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ary  conditions  that  exhibit  strong  velocity  distortion  as  well  as  asymmetry  in  the  turbulence  char¬ 
acteristics.  A  large  amount  of  experimental  data  is  currently  being  processed  to  examine  the 
frequency  dependence  of  the  modes,  the  modal  shape  functions  as  determined  from  the  eigenvec¬ 
tors  of  the  two  point  sturctural  correlation  matrix,  and  finnally  near- tie  id  and  far- fie  Id  sound  radi¬ 
ation  that  is  generated  by  the  shell  vibrations.  These  data  will  be  compared  to  ongoing  analytical 
work  at  Notre  Dame  that  is  focused  on  the  swirling  wake  distribution  downstream  of  the  rotor  and 
the  interaction  of  the  shell  modes  with  the  convecting  pressures. 

In  closing,  the  project  is  complete  and  has  been  successrul  in  meeting  the  original  objectives  of 
creating  the  experimental  tools  and  expertise  that  allows  practical  structural  acoustic  problems  to 
be  studied  using  air-anaiog  experiments.  The  processing  and  presentation  of  the  complete  data  set 
including  acoustic  forcing,  simple  hydrodynamic  forcing,  and  ducted  rotor  forcing  is  still  in 
progress  as  part  of  the  ongoing  ONR  funded  project  involving  inlet  distortion  and  wake  develop¬ 
ment. 
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